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ABSTRACT 


A  computer  program  for  simulating  submarine  motion  in 
six  degrees  of  freedom  is  developed.  The  simulated 
submarine  is  given  a  capability  of  shifting  ballast.  An 
automatic  pitch  and  depth  control  is  designed  for  the 
submarine  simulated  using  optimal  control  theory.  With  the 
depth  and  pitch  in  automatic  control  a  trim  error  signal  is 
developed  by  comparing  the  parameters  of  the  plant  to  that 
of  a  linearized  model.  This  error  is  used  to  implement  an 
automatic  trim  control  that  shifts  ballast  within  the 
submarine  to  achieve  a  neutral  trim. 
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I-   SUBMARINE  CONTROLS 

\  The  submarine  control  problem  differs  from  that  of  the 
surface  ship  in  that  the  submarine  translation  in  three  axes 
and  rotation  about  each  of  these  axes  must  be  considered.  A 
typical  submarine  implements  control  about  these  axes  with 
the  use  of  movable  plane  control  surfaces,  shifting  of 
ballast  water  and  the  propulsion  system. 


The  submarine  operates  and  is  controlled  in  all  six 
degrees  of  freedom.  To  maneuver,  three  sets  of  plane 
surfaces,  the  propulsion  system,  and  a  means  of  shifting 
ballast  are  used.)  This  thesis  considered  only  control  of 
the  plane  surfaces  and  the  shifting  of  ballast  water.  In 
addition  ,  of  the  three  plane  surfaces  available  only  the 
fairwater  planes  and  stern  planes  were  automatically 
controlled,  i.e.  no  automatic  rudder.  The  three  modes  of 
control  were  depth,  p:j.tch  and  ballast  control.  Course 
control  was  not  considered.  ) 

Depth  control  consists  of  the  sensing  of  any  error 
between  ordered  depth  and  pitch  and  the  actual  depth  and 
pitch.  The  planes  were  used  to  implement  the  control  and 
bring  the  submarine  to  the  ordered  depth  and  pitch  thereby 
reducing  the  error  to  zero  or  near  zero. 

Ballast  control  is  the  sensing  of  other  than  neutral 
trim  of  the  submarine.  Because  of  surrounding  water 
temperature  changes,  taking  sea  water  in  the  ship  or  pumping 
it  out,  the  submarine  may  not  weign  the  same  as  the  volume 
of  water  it  displaces.  Also  too  much  water  forward  or  aft 
in  the  ship  initiates  a  turning  moment  around  the  Y  axis 
that  will  eventually  have  to  be  corrected  by  the  planes. 
This  condition  must  first  be  detected  so  a  corrective  action 
can  be  applied  by  shifting  ballast  within  variable  ballast 
tanks  located  inside  the  submarine. 


The  purpose  of  controlling  depth  is  obvious.  pitch 
control  is  required  largely  because  it  is  coupled  to  depth 
control  and  it  is  sometimes  desired  to  maintain  pitch  at 
other  '  tnan  a  zero  angle.  A  means  of  ballast  control  is 
required  to  surface  the  ship,  maintain  depth  under  certain 
conditions,  allow  deptn  control  to  be  accomplished  within 
the  limits  of  the  planes'  capabilities  and  to  reduce  drag 
caused  by  a  steady  state  non  zero  planes  angle.  The 
corrective  action  for  the  last  two  out  of  trim  conditions  is 
normally  accomplished  by  use  of  a  pump  to  shift  water  into 
or  out  of  variable  ballast  tanks  in  the  submarine. 

The  method  of  achieving  an  automatic  control  was  to 
first  develop  an  automatic  controller  that  maintained  an 
ordered  depth  and  pitch.  The  controller  was  to  counteract 
external  forces  by  movement  of  the  fairwater  and  stern 
planes.  These  forces  can  be  caused  by  an  out  of  trim 
condition  or  by  the  ship  dynamics  which  at  different  speeds 
would  cause  the  submarine  to  cnange  depth  in  varying  degrees 
plus  sucn  things  as  sea  state,  etc.  The  controller  designed 
is  unique  to  the  submarine  considered  in  ref (1)  but  the 
method  is  general  and  may  be  applied  to  most  situations. 
With  the  depth  controller  in  operation  the  automatic  control 
of  ballast  was  then  solved.  Again  the  solution  is  unique  to 
the  submarine  considered.  Finally  with  all  automatic 
controls  in  operation  various  runs  were  made  to  determine 
the  capability  of  the  system. 

However,  before  a  controller  design  can  be  started  a 
computer  model  of  the  submarine  must  be  developed. 
Reference  1  contains  a  FORTRAN  program  used  by  NSRDC.  A 
program  language  that  was  more  easily  manipulated  in  terms 
of  filters,  integration  routines,  etc.  was  DSL.  To  use  DSL 
the  equations  of  motion  would  have  to  be  solved. 
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II.   MODELING  A  SUBMARINE  IN  A  DIGITAL  COMPUTER 

A.    EQUATIONS  OF  MOTION  FOR  A  SUBMARINE 

The   equations   of  motion  in  six  degrees  of  freedom  are 
derived  by  equating  the  summation  of  the  forces  and   moments 
to  zero. 
*Z  Forces=0 
Z  Moments=0 

Where  the  vector  force 
F  =  d/dt  Momentum 

=d/dt  (mVe)  =  md/dt  Ve 
m  is  the  mass  and  Ve  is  velocity 
Ve  =  iU+jV+kW 
The  direction  of  positive  motion  is  indicated  in  fig.  II- 1. 

Moment  =  d/dt  (Angular  momentum)  at  the  center  of  gravity. 
Moment  =  d/dt  (iPIxx+ jQIyy  +  kRlzz) 

Where  Ixx  is  the  moment  of  inertia  about   the   X   axis   etc. 
The  direction  of  positive  motion  is  shown  in  Fig.  II-2. 
A  simplified  derivation  cf  the  equations  of  motion  is   found 
in  ref.  2. 

The  equations  are  referred  to  the  body,  i.e.  the  origin 
of  the  coordinate  system  is  located  in  the  body  and  at  the 
center  of  gravity.  Figure  II-3  is  a  submarine  placed  on  the 
coordinate  axis  with  arrows  pointing  in  the  direction  of 
positive  motion.  In  addition  the  movement  of  the  rudder, 
stern  planes  and  fairwater  planes  are  also  depicted  on  Fig. 
II-3  with  the  arrows  shown  pointing  in  the  direction  of 
positive  motion.  The  specific  equations  of  motion  for  a 
submarine  were  developed  for  NSRDC  and  presented  in  ref.  3. 
If   the   coordinate   origin   is  put  at  the  center  of  gravity 
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these  equations  become  the  the  equations   of   ref.   1.    The 

equations   of   motion   for   a   submarine   in   six  deqrees  of 

freedom  are  repeated  from  ref.  1  in  the  Appendix.    For   use 

in    the    computer    programs   that   follow   the   equations 

containing  I  ,1   and  I    had   both   sides   of   the   equation 
x   y       z 

5 

divided   by   1  „    Where   1   is   the   ship   length.    In  the 

3 

equations  containing  m,   both   sides   were   divided   by  \.P 
was  assigned  the  value  of  2. 


The  equations  contain  forcing  terras  for  the  planes  and 
rudder,  shifting  of  ballast  water  and  the  propeller.  The 
propeller  input  is  a  function  fitted  to  a  curve.  It 
contains  the  command  speed  and  actual  speed  as  ratios. 
Normally  there  are  different  sets  of  coefficients  for 
different  propeller  modes  such  as  backing,  slowing  down  etc. 
The  only  mode  considered  was  forward  propulsion  at  a 
constant  speed.  JIt___was___noted  that  the  actual  speed-w-a-s- 
always  slightly  lower  than  the  ordered  speed.  Uc  (command 
speed)  and  U  (actual  speed)  are  in  feet  per  second. 


All  angles   are   in   radians.   Mass   is   in   slugs   and 
distance  in  feet. 


B. 


COMPUTED  SOLUTION  TO  THE  EQUATION'S  OF  MOTION 


To  solve  the  resultant  differential  equations  ref.  1 
partitions  the  equations  into  different  sums.  The  equation 
for  U  contains  no  other  acceleration  term  so  it  is  solved 
separately.  V,P  and  B  equations  all  contain  mutual 
accelerations  so  they  are  solved  using  Cramers  Uule  for  a  3 
x  3  system.  W  and  Q  contain  mutual  accelerations  so  they 
are  grouped  and  solved  as  a  2  x  2  system. 
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The  approach  chosen  in  this  thesis  was  to  solve  the 
whole  system  of  equations  at  once  using  Cramers  Rule  for  a  6 
x  6  system.  This  was  more  methodical  and  the  determinants 
were  solved  only  once  in  the  beginning  of  the  program  using 
the  "INITIAL"  capability  of  DSL. 

Tne  general  forjn  of  the  equations  is,  using  two  degrees 
of  freedom  for  an  example 


K     K 
it    e  1 


K     K 

12     ZZ 


X  =  F  - 


M   S+N      M   S  +  N 

11      11       21     21 


H   S+N     a      S+N 

12      12     22     22 


To  simplify  the  manipulation  let 


I  =  F  - 


M   S+N      M   S+N 

11     11       21     21 


H   S+N      a       S+N 
12     12      zz  zz 


X 


then  by  substituting  with  I 


K  K 

11  2  1 


K  K 

12  22 


X    =    I 


r>J  /V 


defining 
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A  = 


K      K 

11       2  1 


K      K 

12       22 


Using  Cramers  rule  to  solve  for  X   and  X   , 

3  12 


I     K   S 

1       2  1 


I     K   S 

2       ZZ 


toting   the   cofactor  of  K    =  K    and  the  cofactor  of  K 

1122  1 


-K   etc 

21 


X   =1   cofactor  K    +  I    cofactor  K 

11  112  12 


Z 

S  A 


X   =1   cofactor  K    +1   cofactor  K 

Z  Z  ZZ  X  2  1 

2 

s  A 


m 


The  same  general  approach  can  be  used  for  higher  order 
systems  with  the  number  of  cofactors  taken  equal  to  the 
number  of  variables  to  be  solved.  In  the  equations  of 
motion  I  will  be  a  number  of  time  varying  nonlinear 
functions. 

Reference  4  used  the  library  subroutine  DTERM  in 
modified  form  to  solve  similar  equations  for  the  CSMP 
language.  This  method  was  further  modified  to  work,  with  DSL. 
Six  by  six  determinants  were  solved  by  the  program  to  get  a 
solution  for  the  equations  in  the  form  of  the  previous 
paragraph.  By  substituting  columns  and  rows  of  ones  and 
zeros  in  the  row  and  column  wnere  the  variables  I  should  be 
the  cofactor  was  found.  For  example  in  a  3  by  3  case 


0 
K 


22 


23 


32 


33 


=  K   K     -  K   K    =  cofactor  K 

22   3  3         32   2  3  11 


K 


1  1 


K 


31 


K 


12  32 

0       1     0 


=  K   K     -  K   K    =  cofactor  K 

1231         1132  2  3 


Program  #1  was  used  to  solve  for  the  cofactors  in  this 
manner.  Once  the  cofactors  were  found  subsequent  programs 
had  the  cofactors  read  in  as  constants.  In  the  program  the 
cofactors  are  constants  and  the  variable  I  is  dependent  on 
time  and  the  dynamics  of  the  submarine. 
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C.    SUBMARINE  SIMULATION 

The  computer  used  for  the  simulation  was  the  IBM  360/67 
located  at  the  W.R.  Churcn  Computer  Facility,  Naval 
Postgraduate  School.  The  programming  language  was  the  Naval 
Postgraduate  School  version  of  DSL/360.  Computer  program  #1 
is  a  simulation  of  a  submarine  in  motion  at  15  knots.  The 
coefficient  values  used  are  obtained  in  ref.  1.  As  was 
previously  discussed  the  equations  of  motion  are  derived  in 
the  body  fixed  system.  A  trigonometric  conversion  of  the 
velocities  is  used  in  the  auxiliary  equations  (appendix)  to 
transform  the  equations  to  the  eartn  reference  system  with 
the  origin  at  x(to),  y  (to)  and  z(to). 

Reference  1  contains  several  maneuvers  used  to  check 
out  the  equations  and  coefficients.  One  maneuver  is  the 
vertical    plane   overshoot   where   the   stern   planes   are 

o  o 

deflected  to  +20  lf    when  the  ship   is   pitched   to   '5     the 

stern   plane  angle  is  reversed.  Program  #1  executes  the  same 

maneuver  with  fig.  II-4  and   II-5   showing   depth, pitch   and 

stern  plane  angle.  The  results  compared  favorably  with  those 
of  ref.  1. 

The  original  equations  contained  forcing  terms  that 
were  to  represent  water  being  blown  from  the  main  ballast 
tanks.  In  the  programs  developed  for  this  thesis  tue 
forcing  terms  represented  water  being  moved  in  variable 
ballast  tanks  with  the  capability  of  shifting  ballast 
between  tanks  or  between  any  tank  and  sea  . 

To  control  the  trim  of  a  submarine  a  means  of  shifting 
ballast  water  must  be  provided.  For  the  computer  simulation 
a  trim  pump  must  also  be  simulated.  Also  trim  tanks  must  be 
positioned  in  the  ship  model.  The  position  of  the  tanks  was 
an   arbitrary  choice    with  some  traditional  basis.   The  tanks 
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were  positioned  as  indicated  in  fig.  I3>6. 

Many  options  are  available  in  the  shifting  of  variable 
ballast  water.  The  options  depend  on  the  number  of  tanks 
and  their  individual  purpose.  To  simplify  the  computer 
simulation  the  options  will  be  decreased  considerably .  Only 
three  tanks  will  be  considered.  The  auxiliary  tank  is  placed 
at  the  center  of  gravity  and  will  be  used  to  adjust  any 
overall  mass  error.  That  is  a  weight  differential  between 
the  snip  and  the  water  it  displaces.  The  forward  and  after 
trim  tanks  will  be  used  only  to  adjust  ballast  within  the 
ship.  The  result  will  allow  negative  values  for  the  content 
of  a  tank.   Tne  options  are, 


1.  Heavy-  pump  auxiliary  tank  to  sea. 

2.  Light-  flood  auxiliary  tank  from  sea. 

3.  Heavy  forward-  pump  from  forward  trim  tank  to  after 
trim  tank. 

4.  Light  forward-  pump  from  after  trim  tank  to  forward 
trim  tank. 

Any  ballast  correction  will  be  accomplished  by  first 
correcting  the  overall  mass  error  then  achieving  a  forward 
and  aft  trim. 

Using  the  pump  from  Program  #2  with  Program  #1  the 
effects  of  shifting  ballast  was  demonstrated  in  the 
following  simulations*  The  planes  are  put  at  zero  and  not 
moved  during  the  run.  Figure  II-7  is  the  resultant  pitch 
change   when   2C/00   pounds   of   ballast   is  flooded  into  the 

auxiliary  tank  at  a  rate  or  3860  lbs./min.   Because   of   the 

o 

ship   dynamics   a   2.7    up   angle   is  taken  even  though  the 

auxiliary  tank  is  located  at  the  center   of   gravity.    With 
the   planes  left  at  zero  the  pitch  angle  remains  constant  at 
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this  angle.  Figure  II-8 ) produces  a  surprising  result,  with 
2000  pounds  flooded'  in  the  ship  went  up  240  feet  in  3.5 
minutes.  (A  negative  number  for  depth  indicates  up  and  a 
positive  number  indicates  down.)  However  keeping  in  mind 
that  the  submarine  assumed  an  up  angle  the  resultant  depth 
change  makes  more  sense.  Figure  II-9  was  plotted  with 
ballast  pumped  out  of   auxiliary   tank   to   sea.    The   ship 

o 

assumes   a   .9    degree  up  angle.   The  depth  is  seen  in  Fig. 

11-10  to  change  by  135  feet  in  3.5  minutes.  Figure  11-11 
helps  to  resolve  the  apparent  problem.  In  fig.  11-11  no 
ballast  was  transferred  yet  the  ship  assumed  an  up  angle  of 
significant  degree  with  the  planes  on  zero.  Reference  5 
explains  that  most  military  submarines  will  generate  a 
positive  pitch  due  to  the  ship  structure  if  the  planes  are 
zeroed.  The  submarine  in  ref.  1  is  fictitious  and  the 
coefficients  were  randomly  selected  but  it  still  appears  to 
have  some  of  the  expected  characteristics.  Finally  fig. 
11-12  and  13  were  generated  with  ballast  shifted  from  after 
trim  to  forward  trim. 


The  model  performed  as  expected  and  the  trim  pump  and 
planes  produced  the  desired  result.  The  next  problem  solved 
was  the  designing  of  automatic  controls. 
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Figure  II-l 
Axes  Definitions 


Figure  1 1 -2 
Axes  Definitions 
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Pitch  vs.    Time 
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Pitch  vs.   Time 
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III.  AUTOMATIC  PITCH  AND  DEPTH  CONTROL 

A.    OPTIMAL  SOLUTION  TO  THE  LINEARIZED  SUBMARINE 
EQUATIONS 

Automatic  depth  and  pitch  control  can  be  achieved  by 
various  methods.  For  example,  because  of  the  fairwater 
planes  having  minimal  effect  on  pitch  a  frequently  used 
combination  is  to  control  the  depth  witn  the  fairwater 
planes  and  pitch  with  the  stern  planes.  Another  way  is  to 
set  a  fixed  ratio  between  stern  planes  angle  and  fairwater 
planes  angle.  In  this  case  pitch  is  not  "controlled"  but  is 
used  to  achieve  ordered  depth  and  the  planes  are  used  to 
change  pitch.  The  approach  selected  in  this  thesis  was  to 
let  the  planes  individually  achieve  what  ever  angle  was 
reguired  to  meet  an  optimal  control  based  on  a  cost  function 
to  be  described  later. 

Because  of  the  nonlinearity  of  the  equations  of  motion 
a  linearized  model  was  reguired  to  solve  the  equations 
required  for  an  optimal  control  system.  The  required  gains 
were  then  determined  using  the  linearized  equations. 
Linearizing  the  equations  about  some  operating  point  was 
found  to  be  impractical  because  the  operating  point  can  not 
be  constant  and  the  nonlinearities  in  the  equations  present 
a  formidable  problem  in  taking  derivatives.  Linearization 
was  accomplished  by  dropping  all  nonlinear  terms  with  the 
following  justifications.  Referring  to  appendix  A  all  terms 
involving  Hi  are  dropped  because  the  linearized  model  is  in 
trim.  N'=1  because  Uc=U.  U  is  a  measured  quantity  and  the 
only  terras  of  interest  are  H  and  Q,  i.e.  the  problem  is 
reduced  to  two  degrees  of  freedom.  The  last  result  is 
obtained  because  only  the  depth  and  pitch  are  parameters  to 
be  controlled.  Simple  truncation  of  the  equations  of  motion 
will  eliminate  the  other  nonlinear   terms.    The   linearized 
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equations  of  motion  are  as  follows 


m-Zw   -lZg   i  ■  W 


-Mw/1    Iy-Mg 


U 


UZw/1       OfZqj  1  "  W 


UMw/i 


UMg/1    Q 


2  2 

U  Zds/1     U  Zdb 


2         Z  Z 

U  Mds/1     U  Mdb/1 


Ds 


Db 


let 


m-Zw    IZq 


Htf/1   Iy-Mq 


-3 


=Inv 


=  Inv 


UZw/1     UZq 


UMw/1     UMq/1 


+  Inv 


U  Zds/1     U  Zdb/1 


2         2 

U  Mds/1 
u. 


2        2 

U  Mdb/1 


Ds 


Db| 
U  -J 


The  values  of  the  inverted  matrix  are  constants  and 
when  solved  for  can  be  directly  substituted  as  constants. 
The  equations  were  then  written  in  a  acre  general  fora  that 
would  lend  itself  to  coding  for  a  digital  computer  program 
as  follows. 


A     A 

11      2  1 


A     A 

12      22 


3      3 

11       2  1 


B      3 

12       22 


Ds 


Db 
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let 


Determ=[  (Iy-Mg)  (m-Zw)-Mw  Zq  J 


then 


A   =[  (ly-Mq) Zw  +  Zq  WwJU/lDeterm 


A   =£MwZw+ (m-Zw)  Mw  ]0"/l  Determ 


A   =[  (Iy-Mq)  Zq+Zcj  MqJU/Determ 


A   =[  MwZg+ (m~Zw) Hg JU/lDeterm 

2  2 


B   =[  (Iy-Mq) Zds  +  ZqMds]U  /IDeterm 

2    2 

B   =[  MwZds+  (c-Zw)  Mds  JU  /l  Determ 

12 


B   =[ (Iy-Mq) Zdb+Zq  Mdo]U  /IDeterm 

2    2 

B   =[MwZdb  ♦  (m-Zw)  MdbJLJ  /l  Determ 
zz 


The  problem  was  treated  as  a  linear  tracking  problem 
with  the  control  designed  around  the  linearized  model-  The 
methods  of  ref.  6  were  used  for  an  optimal  control  with  tne 
state  equations  in  the  general  form 


X  =  AX  ♦  BU 


which  represents 
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Ds 

• 

Q 

=A 

Q 

♦  B 

Db 

—  - 

* 

_      _ 

X  = 


and 


U  = 


Ds 


Db 


A  and  B  are  as  previously  defined. 


The  cose  function  to  be  minimized  is 


.tf    <j 


J- 1/2   /.  [  E    QE    +    U    3U]dt 


/V  •vz       ^/v/ 


In  the  derivation  of  the  optimal  control  equation  the 
control  returns  the  system  to  the  origin.  In  this 
controller  r  was  not  restricted  to  the  origin  so  X  was 
equivalent  to  E  where 


E  =  X  -  r 

*\S  rss  /\s 


and  r  is  the  ordered  depth  and  pitch.  A, 3,0,1*,  and  r  are  all 
constant,  X  and  U  are  time  varying  and  the  plant  is  assumed 
to  be  completely  controllable*    (The   last   assumption   was 
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supported   i>y   the  solution  results.)  T   is  considered  fixed 

f 

and  X  (t  )  is  free.   The  Hamiltonian  is  given  by 

<~   f 


T       T       T       T 
Ham  =  1/2[E  QE    +  U  RU  ]  +p  AX  *    p  BU 


The  costate  equation  is 


p  =^Hara/^X  =-QX-A  p  ♦  Qr 
=*QE-A  p* 


The   asterisks  indicate   optimal  controls  and  trajectories, 
The  optimal  control  is  then  found  by 


«^Ham/c)U  =  BU  ♦  B  p*  =  0 


r^r*/  rsj        ^j 


*    -1  T  * 
0  =-fi   B  p 


!             -1     T 
A        ■     -BR       B 

-     - 
*• 
E 

!          T 
-Q     .     -A 

P 

From  ref.  5  the  boundary  conditions  arc  p  (t  )  =  0 

/N/       f 

defining  p  =KE 


taking  the  derivative  of  p  and  dropping  the  stars 


p  =KE+KX 


/v«v  ^v  ^ 
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by    substitution 

in 
p    =  KE    ♦    KX    =    -QE-A    p 

*f  fs//\f  rvj  /v/v      /v      pj 


-I  T 
X  =  AE-BE   B  p 


p  =  KE 


then 


-IT  T 

K=-KA  +  KBB   B  K  -  Q  -  A  K 


and 


U  =-R   B  p  =  -B   B  KE 


As   descriced   in   ref.   5  K (t  )  =  0.0.   K  is  symmetric  i.e. 
K   =K    etc.  so   n(n+1)/2   differential   equations   must   be 

12    21 

solved  where  n=4.   The  solution  of  the  differential  equation 
for  K  will  provide  the  optimal  gains  for  U. 


Referring  to  the  original  equations  and  putting  them  in 
state  variable  form 
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R  and  Q  are  the  weighting  matrices. 

The    following    equations    result   from   the   matrix 
manipulation. 
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F  =    B       K  +    B       K 

11  1112  1214 

F  =    B       K  +    B       K 

21  1122  12     2* 

F  =    B       K  +    B       K 

31  1123  1234 

F  =    B       K  +    B       K 

♦1  1124  1244 

F  =    B       K  ♦    B       K 

12  2112  2214 

F  =   B       K  +    B       K 

22  2122  22      24 

F  =    B       K  ♦    B       K 

23  212  3  22      34 

F  =    B       K  +    B       K 

2*  2124  224  4 


2  2 

K         =    -E    +  (F  +    F         )/C 

11  11  12 


K  =    -(K          +    K       A          +  K       A       )  *     (F       F          +F         F)/C 

12  11                   12      11  14      12                         21       11                  22          12 

K  =     (F       F          +    F       F  )  /C 

13  3111  32      12 

K  =    -<K       A         +    K          +  K      A       )  *     <F      F         +FF)/C 

14  12       12                  13  14      22                         41      11                  42      12 


22 


-2(K         +    K       A         +    K      A      )     +     (F  +    F  ) /C 

12  22      11  24       12  2  1  22 


K         =-(K  +KA  +    K       A       )     +     {F       F  +F         F/C 

23  13       23   11       34   12  31   21       32         22 

K    =-(K    +  K_  A    +  K   A   )  -  (K   A    +K    +KA) 


24 


14 


24   11 


44   12 


22   2  1 


23 


24   2  2 


t 


(F   F    +  F   F   )  /C 

4  121        42   22 


K  =  (F     +F    )  /C-D 

3  3         3  1        32 

K  =  -  (K   A    +  K    +  K   A   )  ♦  (F   F    +  F   F   )  /C 

34  23      21  33  34      22  41      31  42      32 

2  2 

K  =    -2(K       A         ♦    K  +    K       A       )     +     (F  ♦    F  ) /C 

44  2421  34  4422  41  42 


Since  only  the  value  of  K (t  )  is  known  these  equations   oust 

~   f 
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be  solved  backwards  i.e.  from  t   to  t  .   To  do  this  let 

f     o 

r  =t  -t 
f 

df=  -dt 

when  t=t   ,T=0 
f 

K  =  dK/dt  =  -dK/dT 


The  signs  of  all  the  equations  are  reversed  and  the   initial 

conditions  are  those  conditions  at  t  .   Up  to  this  point  the 

f 

solution  is   general   and   when   solved   the   solution   will 

provide  gains  for 


-1  T 

U  =  -R      B  KE 


=-1/C 


F  F  F 

11     21      31 


F  F  F 

12     22     32 


4  1 


♦J 
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The  calculations  of  the  gains  that  follow  apply  only  to 
the  submarine  in  ref.  1.  The  submarine  described  in  ref.  1 
is  fictitious  and  used  for  demonstration  purposes  only.  DSL 
again  is  used  to  solve  for  the  gains.  The  program  used  is 
program  #2  which  also  contains  the  controllers  designed  in 
the  following  sections. 

Prograjtt  #2  first  calculates  the  steady  state  gains  used 
in  the  controller.  While  the  gains  are  calculated  the  ship 
simulation  is  held  at  zero.  After  a  specified  time  the 
simulation  starts  with  the  gains  determined  by  the  specified 
weighting  matrix. 

Selection   of   the   weighting   matrix  to  be  used  can  be 
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started  off  with  some  judgment  as  to  the  severity  of 
different  errors.  For  instance  the  pitch  equations  are  in 
radians  and  the  depth  equations  are  in  feet.  It  is  obvious 
that  a  one  radian  pitch  error  is  more  severe  than  a  one  foot 
depth  error.  This  immediately  would  lead  to  a  higher  weight 
on  the  pitch  error  signal.  It  was  found  that  if  too  much 
emphasis  was  placed  on  pitch  the  depth  error  in  steady  state 
would  become  unacceptable.  This  lead  to  a  trial  and  error 
process  to  select  an  acceptable  weighting  matrix  combination 
that  would  give  the  desired  results. 

To  operate  as  an  optimal  controller  the  gains  would 
have  to  be  generated  on  line.  This  would  require  a 
significant  computer  capability.  In  an  infinite  duration 
process  the  gains  would  be  constant  as  demonstrated  in  fig. 
III-1  through  7.  The  optimal  control  can  then  be 
approximated  witn  a  steady  state  gain  if  the  control 
interval  is  long  with  respect  to  the  decay  time  of  the 
optimal  gains.  Figure  III-1  through  III-7  indicate  a 
relatively  short  rise  and  fall  time  for  the  gains  just 
before  t  . 


As  determined  by  observing  the  equations  for  K  the 
gains  are  functions  of  the  velocity  U.  A  number  of 
different  runs  were  made  at  different  velocities  and  the 
following  approximate  proportionalities  were  observed.   K  , 


K   #K    and  K    are  inversely  proportional  to  U  .    K   ,K   , 

14    23  34  _   r  ZZ         24 

3 

and   K    are   inversely   proportional  to  U  .   The  runs  were 

4  4 

made  at  25.33  feet  per  second  so   the   following   adjustment 
was  made  to  the  gains  used  in  the  programs. 
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2         2 

ZM1=(25. 33)   K   /U 

12 
2  2 

ZM2=(25. 33)   K   /U 
l  ♦ 


3  3 

ZM3=(25. 33)   K   /U 

22 


2  2 

ZM4=(25.33)   K   /* 

23 


ZM5=(25. 33)   K  /U 

24 

2  2 

ZM6=(25. 33)   K    /U 

34 

3  3 

ZM7~(25. 33)   K   /U 

4  4 


The   new   gains   are   given  by  ZM1  through  ZM7.  The  B  matrix 

2 

contains  U   so  a  cancellation  of  U  takes  place  in  the   final 
eguation. 

B.    AUTOMATIC  PITCH  AND  DEPTH  CONTROL 

Figure  III-8  and  III-9  were  simulations  with  C=1.0, 
D=100.0  and  E=1.0  or  the  control  input  was  weighted  to  1.0, 
the  pitch  error  was  weighted  to  100.0,  and  depth  error 
weighted  to  1.0.  The  ordered  depth  change  was  10.0  feet.  The 
controller  initiated  the  proper  planes  but  because  the  pitch 
was  so  large  when  ordered  depth  was  achieved  the  overshoot 
pushed  the  control  system  into  instability.  Also  the  ship 
pitches  in  the  wrong  direction.  The  initial  plane  response 
is  to  a  depth  error  and  both  planes  were  positioned  as  seen 
in  fig.  111-10.  The  ship  pitched  up  with  an  angle  too  large 
to  overcome. 

The  next  run  was  made  after  weighting  the  pitch  and  the 
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control  input  more  heavily.  The  weighting  matrix  used  was 
C=10.0,  0=300.0,  and  E=1.0.  Figure  111-11  is  the  resultant 
depth  change  tor  a  10  foot  step  input  in  ordered  depth.  The 
resultant  steady  state  depth  error  was  .29  feet  which  was 
acceptable.    As   before   the   ship   initially  pitched  up  to 

o 
about  12   and  then  returned   to   a   steady   state   angle   of 

o 
+1.28    as   seen   in   fig.   111-12.   Figure  111-13  shows  the 

fairwater  and  stern  plane  angles  used   to   achieve   the   new 

o  o 

ordered   depth.   The  maximum  -80   stern  plane  angle  and  -50 

fairwater  plane  angle  are  unacnievable  on  a   real   ship   and 

the   rates   of  change  of  depth  and  pitch  were  unnaccep table. 

When  a  limit  constraint  was  placed  on  the  plane   angles   the 

system   went   into  a  limit  cycle.   The  steady  state  error  in 

pitch  was  large  enough  to  prevent  the  txim  controller  in  the 

next   section   from  detecting  a  ballast  change  to  any  degree 

o 
of   accuracy.    Lastly   the   pitch   variation   of   12     was 

considered  unacceptable  for  only  a  10  foot  depth  change. 

Larger  emphasis  was  required  on  pitch  but  this  caused 
the  plane  angles  to  become  larger  and  if  the  weighting  on 
plane  angles  was  increased  the  system  became  insensitive  to 
either  depth  error  or  pitch  error  depending  on  the  weight. 
The  problem  appeared  to  be  a  result  of  the  initial  depth 
error.  To  limit  this  effect  with  a  large  step  input  the 
depth  and  pitch  error  were  limited  as  shown  in  fig.  111-14. 
Also  a  larger  emphasis  was  required  on  pitch  error  in  the 
weighting  matrix.  The  weighting  matrix  chosen  was  C=10.0, 
D=3000.0#  and  E=1. 0. 

In  fig.  111-13  a  significant  noise  level  is  seen  in 
the  plane  positions.  To  reduce  this  noise  the  plane  angle 
ordered  was  filtered  with  a  simple  low  pass  filter. 

The   tests  were  run  again  and  this  time  a  20  foot  depth 
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change  was  ordered.  The  magnitude  of  the  saturation  limit 
can  be  adjusted  for  different  results  .   The  values  selected 

o 

for  saturation  were  ±  two  feet  for  depth  error  and  ±  10   for 

pitch   error.   Major  differences  to  be  noted  in  the  runs  are 

the  reduced  rate  of  depth  change.   A  20   foot   depth   change 

o 
took   about   45  seconds.   The  maximum  pitch  assumed  was  -.8 

The  steady  state  depth  error  was  .2  feet  and   the   steady 

o 
state   pitch   error   was   . 14  .   This  run  was  plotted  in  fig 

111-15,  16  and  17. 

The  controller  was  tested  for  its  ability  to  maintain 
depth.  To  check,  this  the  ordered  depth  was  put  at  zero. 
The   controller  maintained  a  steady  state  depth  error  of  .16 

o 
feet  and  a  steady  state  pitch  error  of  .14    as  seen  in  fig. 

111-18   and   111-19.   The  plane  angles  are  presented  in  fig. 

111-20.  Note  that  a  steady  state  fairwater   plane   angle   of 

o 

-4.0    was  reguired  to  maintain  depth.   What  this  implies  is 

that  to  run  the  simulated  submarine  at  a  zero  pitch  and  on 
the  ordered  depth  with  zero  plane  angles,  ballast  would  have 
to  be  added. 

The   depth   controller   was   also   reguired  to  maintain 

o 

depth  in  a  turn.   A  simulated   rudder   angle  of   10      was 

added.  The  rudder  motion  is  not  restrained  in  rate  for  this 
simulation.  After  50  seconds  the  rudder  angle  was  put  to 
zero.  Figure  111-21  plots  the  roll  developed  as  a  result  of 
the  manipulation.  The  transient  in  pitch  and  depth  was  small 
as  indicated  in  fig.  111-22  and  23.  In  particular  fig. 
111-24  shows  a  positive  fairwater  and  stern  plane  angle 
reguired  as  the  submarine  settles  in  the  stern. 

In  early  runs  it  was  apparent  that  the  submarine  was 
"light"   and   ballast   should  be  added.   This  was  based  on  a 
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0 

steady  state  fairwater  plane  angle  of  about   -4.0     with   a 

near   zero   depth   and   pitch   error.    In  the  following  run 

ballast  was  added  to  the  auxiliary  tank  at  about  8000  pounds 

per   minute  for  2  minutes  while  the  planes  were  in  automatic 

control.   Figure  111-^25  and  26  show   a   steady  state   depth 

o 
error   and   pitch   error   of   about   -.02   feet  and  .04   as 

o 

compared  to  .16  feet  and  .14    with  out  the   ballast.    Also 

the  steady  state  fairwater  plane  angle,  Fig.   111-27  ,  shows 

o  o 

a  steady  state  angle  of   about  -1.5    compared   to   -4.0  . 

This   confirms   that  ballast  should  be  added  to  run  the  ship 

in  •■trim". 


The  depth  controller  was  operating  satisfactorily  with 
a  controller  that  is  somewhat  less  than  optimal  for  large 
transients  and  near  optimal  for  small  transients.  The  last 
problem  to  solve  was  to  determine  ballast  error  in  the 
submarine  and  correct  it. 
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IV.  AUTOMATIC  TRIM  CONTROL 

A.    METHODS   OF  TRIM  CONTROL 

A  continuous  monitoring  and  control  of  the  submarine 
trim  is  required.  A  submarine  that  is  "heavy"  runs  with  the 
fairwater  planes  on  a  constant  positive  angle.  If  the 
submarine  is  too  heavy  the  fairwater  planes  angle  reaches 
saturation  and  a  large  positive  pitch  is  required  to 
maintain  depth.  The  speed  of  the  submarine  along  with  the 
angle  of  inclination  of  the  planes  determine  the  lifting 
force  developed  by  the  planes.  At  a  high  speed  a  submarine 
that  is  "heavy"  may  be  well  within  the  control  limits  of  the 
planes.  When  the  speed  is  decreased  the  decreased  lifting 
forces  would  cause  larger  plane  angles  to  be  ordered  until  a 
possible  saturation  of  the  planes  could  occur.  A  submarine 
at  any  speed  with  a  constant  steady  state  fairwater  or  stern 
planes  angle  is  experiencing  a  drag  that  is  reducing  the 
efficiency  of  the  propulsion  plant.  These  conditions  must  be 
corrected  to  efficiently  operate  the  submarine. 

Two  options  for  automatic  trim  control  were  considered. 
One  was  to  determine  the  steady  state  lifting  forces  caused 
by  excessive  ballast  and  display  the  magnitude  of  the 
ballast  error  allowing  a  manual  correction.  The  next  was  an 
automatic  controller  in  the  loop  pumping  or  moving  ballast 
as  needed  to  maintain  trim. 

The  first  option  would  be  accomplished  by  solving  the 
linearized  equations  developed  in  previous  sections  in  the 
form 


Force=AW  +  EQ  ♦  CDs  +  EDb 
Moment=AW  +  EQ  +    CDs  +  EDb 
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Where  in  the  steady  state,  if  the  plane  angles  are  non  zero 
the  forces  and  moments  are  equivalent  to  the  weight  of 
required  ballast  to  be  moved.   This  option  was  not  explored. 

The  second  option  was  accomplished  using  a  model 
following  technique.  The  velocities  of  a  linearized  model 
running  in  a  computer  and  those  of  the  submarine  are 
compared  as  in  fig.  IV- 1,  Both  the  model  and  the  plant  have 
the  same  plane  angle  inputs.  For  zero  plane  angles  and 
operation  on  ordered  depth  tne  velocities  in  the  vertical 
direction  and  around  the  Y  axis  for  both  the  plant  and  model 
are  zero.  Perturbations  caused  by  the  planes  to  achieve  a 
new  ordered  depth  should  cause  both  the  plant  and  the  model 
to  accelerate  in  the  same  direction.  However  if  a  steady 
state  plane  angle  is  required  to  maintain  depth  or  pitch  in 
the  submarine  these  plane  angles  represent  accelerations  in 
the  model.  A  comparison  of  the  submarine  which  is  at  nearly 
zero  velocity  and  the  model  wnich  is  accelerating  would 
indicate  a  ballast  differential  between  a  model  in  "trim" 
and  a  submarine  out  of  "trim".  There  will  be  some  error  in 
the  model  but  for  small  perturbations  the  error  should  be 
small.  The  model  actually  used  was  the  same  as  used  in 
previous  sections  for  the  derivation  of  the  depth 
controller. 

The  need  for  the  watch  officers  to  maintain  absolute 
control  of  the  ballast  shifted  in  or  out  of  the  ship  would 
prevent  an  automatic  trim  controller  from  working  in  other 
than  a  passive  override  mode.  In  this  mode  the  controller 
signals  the  operator  that  corrective  action  must  be  taken  to 
maintain  the  proper  trim.  The  operator,  if  conditions 
permit  will  actively  allow  the  controller  to  perform  its 
task.  With  no  action  on  the  part  of  the  operator  the  system 
would  not  operate. 
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fl.    MODEL  COMPARISON  FOR  TRIM  CONTROL 

In  all  the  simulations  that  follow  the  velocities  used 
are  W  and  Q.  The  error  resulting  is  small  as  long  as  pitch 
and  roll  are  small.  It  was  observed  that  with  automatic 
control  the  pitch  was  kept  relatively  small.  The  effect  of 
the  roll  will  be  discussed  later.  All  runs  are  made  with 
the  submarine  in  an  automatic  depth  control  mode. 

To  get  an  idea  of  the  magnitude  of  the  velocities 
involved  the  first  runs  consisted  of  comparing  the  model  and 
submarine  velocities  in  the  W  and  Q  direction  while  manually 
initiating  a  shift  of  ballast.  In  fig.  IV-2  ballast  was 
pumped  into  the  auxiliary  tank.  The  bottom  curve  is  the 
fairwater  planes  angle  which  moved  towards  a  zero  position. 
This  is  the  same  result  observed  in  the  previous  section. 
The  top  curve  is  the  stern  planes  angle  and  the  middle  curve 
is  10  rimes  the  velocity  error.  The  error  signals  are 
multiplied  to  fit  them  on  the  same  plot  as  the  planes.  Note 
that  the  velocity  error  crossed  zero  at  approximately  the 
same  time  as  the  fairwater  planes  were  at  zero.  This  was 
the  desired  response  to  adjust  trim  for  a  zero  fairwater 
plane  angle. 

This  simulation  run  was  repeated,  this  time  ballast  was 
shifted  from  after  trim  tank  to  forward  trim  tank.  In  fig. 
IV-3  the  fairwater  planes  angle  and  stern  planes  angle, 
curve  1  and  2,  were  adjusted  to  overcome  the  turning  moment 
caused  by  the  ballast  shift.  Curve  3  is  2000  times  the 
velocity  errcr  in  the  Q  direction.  This  curve  crossed  zero 
at  about  the  same  time  as  the  stern  plane  angle  crossed 
zero.  This  was  the  desired  response  for  a  zero  stern  plane 
angle  when  the  sunmarine  was  in  "trim". 

The  controller  used  was  an  on-off  controller.  The  pump 
operated  at  a  fixed  rate.   The  error  detectors  were  to  sense 
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an  error,  determine  the  direction  of  the  error  and  then 
operate  the  pump.  To  prevent  the  pump  from  operating 
continuously  a  deadzone  was  inserted  in  the  error  signal 
channel.  the  deadzone  was  used  to  limit  cycling  of  the  pump 
during  transients  but  its  magnitude  was  constrained  by  the 
plane  angles  allowed  in  steady  state.  Adjustment  of  tne 
deadzone  would  also  determine  the  frequency  of  pump  cycling 
and  the  steady  state  accuracy. 

In  fig.  IV-4  the  vertical  velocity  error  and  the  steady 

state   accuracy   were  plotted  with  a  history  of  maneuvers  as 

o 
follows.   At  time  zero  the  rudder  was  put  at  a  left  10   ,  in 

50   seconds   it   was  moved  back  to  zero.   At  100  seconds  the 

depth  ordered  was  -20.0  feet.   The  peak  observed  at  time   70 

seconds   would  prevent   any  fine  control  of  trim.  This  peak 

was  the  result  of  the  settling  of  the  submarine  in  a  turn  as 

observed  in  previous  sections.   The  settling  or  squatting  is 

compensated  for  with  plane   angle   and   these   plane   angles 

accelerate   the  model.   This  is  not  a  serious  poblem  and  can 

be  solved  in  a  number  of  ways.   One  way  would  be  to  increase 

the   degrees   of  freedom  in  the  model.   This  would  of  course 

be  a  more   difficult   model   to  construct   in  a   computer. 

Alternately  the  system  can  disengaged  in  a  turn.   This  would 

reduce  the  capability  of   the   system   by   a   small   degree. 

Since    the    time    spent   in   actually   turning   is   small 

disengaging  might  be  the  best  course  of  action.   For  use   in 

this  thesis  it  was  assumed  that  this  second  course  of  action 

was  taken.    The   deadzones   were   then   selected   at   .0008 

radians/second   for   rotational   velocity   error   and   .08 

feet/second  fcr  verical  velocity  error. 

The  controller  had  two  separate  modes.  One  was  to  move 
ballast  in  and  out  of  auxiliary  tank  to  sea.  The  other  was 
to  move  nallast  between  forward  and  after  trim  tanks.  These 
modes   were   first   checked   out   separately.     First    the 
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ballast  the  pump  was  shifted  to  pump  from  auxiliary  for  79 
seconds.  At  this  point  the  controller  stopped  pumping  out 
of  the  auxiliary  tank  and  again  pumped  from  after  trim  to 
forward  trim  for  20  seconds  and  then  the  pump  stopped.  70 
seconds  after  the  pump  stopped  it  started  and  pumped  for  20 
seconds  out  of  the  auxiliary  tank.   The  end   result   on   the 

o        o 
steady   state  angle  of  the  fairwater  planes  was  .6    and  .2 

on  the  stern  planes.   Figure  IV-9  shows  the  tank   levels  as 

ballast   is   moved   to  achieve  a  neutral  trim.   Figure  IV-10 

and  11  plot  the  planes  angles.   Tne   planes   move   toward  a 

zero  steady  state  angle  and  a  neutral  trim  is  approached. 

To  denomstrate  the  effect  of  the  deadzone  the   vertical 

velocity   error   deadzone  was  reduce  to  .04  feet/second  from 

.08  feet/second.  Figure  IV- 12  is  the  tank  levels   where   the 

initial  out  of  trim  condition  was  the  same  as  for  fig.  IV-9. 

A  comparison  of  fig.  IV-10  and  11  with   fig.   IV-1J   and   14 

shotfs   the   steady   state   planes  angles  were  closer  to  zero 

with  the  smaller  deadzone.   Figure  IV- 15   and   16   show   the 

depth    controller   while   ballast   is   shifted.    Pitch   is 

o 
controlled  to  .02   error  and  depth  to  a  . 1  foot  error. 

When  the  ship  was  ordered  to  change  depth  by  20  feet 
with  the  ship  in  trim  the  model  followed  closely  enough  so 
that  no  pump  commands  were  given. 

An  additional  feature  in  the  system  was  tested.  On 
some  maneuvers  it  is  desirable  to  operate  the  submarine 
••heavy".  This  is  normally  accomplished  by  simply  adding  for 
example  5000  lbs  to  the  auxiliary  tank.  Whether  or  not  the 
submarine  is  actually  5000  lbs  heavy  depends  on  the  initial 
state  of  the  trim.  If  the  submarine  is  initially  2500  Ins 
light  it  is  now  only  2500  lbs  heavy  vice  the  desired  5000 
lbs.  By  adding  a  weight  forcing  term  to  the  model  in  the  W 
equation  and  assigning  a  value  to  this  forcing  term  when  the 
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heavy  or  light  conditon  is  desired  the  pump  will  shift  water 
into  or  out  of  the  auxiliary  tank  until  the  desired 
condition  is  reached.  The  accuracy  of  accomplishing  this  is 
determined  by  the  deadzone  in  the  vertical  velocity  error 
signal.  For  the  submarine  simulated  a  value  of  .02 
feet/second  worked  nicely. 

The  program  used  for  the  simulation  is  program  #2.  Once 
the  gains  and  cofactors  are  determined  they  can  be  read  into 
the  program  as  parameters  and  the  subroutines  and  Riccati 
Equations  (gain  equations  on  page  38)  can  be  removed.  The 
controllers  with  filters  as  used  is  shown  in  fig.  IV-17. 
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Model  Comparison  Block  Diagram 
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Figure  IV-1 
Model  Comparison  Block  Diagram 
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Figure    IV-2 
Curve    1,    Fairwater   Planes  Angle   vs.    Time 
Curve   2,    Stern   Planes  Angle   vs.    Time 
Curve    3,    Vertical  Velocity  Error  vs.    Time 
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Figure  IV-3 
Curve  1,  Fairwater  Planes  Angle  vs.  Time 
Curve  2,  Stern  Planes  Angle  vs.  Time 
Curve  3»  Rotational  Velocity  Error  vs.  Time 


81 


o 


o 

<D 

W 

\ 

-P 
Q) 
(D 

CO 
O 


-0-00 


zs.oo 


(X10  ) 

50  Seconds/inch 


Figure  IV-4 
Vertical  Velocity  Error  vs.  Time 
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Auxiliary  Tank  Level  vs.    Time 
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Figure  IV-6 
Curve  1,  Fairwater  Planes  Angle  vs.  Time 
Curve  2,  Stern  Planes  Angle  vs.  Time 
Curve  3,  Vertical  Velocity  Error  vs.  Time 
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Forward  Trim  Tank  vs.    Time 
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Figure   IV-8 
Curve   1,    Fairwater   Planes  Angle  vs.    Time 
Curve   2,    Stern   Planes  Angle    vs.    Time 
Curve    3,    Rotation  Velocity  Error  vs.    Time 
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Curve  1,  Auxiliary  Tank  Level  vs.  Time 
Curve  2,  Forward  Trim  Tank  Level  vs.  Time 
Curve  3,  After  Trim  Tank  Level  vs.  Time 
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Figure    IV-10 
Fairwater  Planes  Angle   vs.    Time 
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Figure   IV-11 
Stern   Planes  Angle  vs.    Time 
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Figure  IV-12 
Curve  1,  Auxiliary  Tank  Level  vs.  Time 
Curve  2,  Forward  Trim  Tank  Level  vs.  Time 
Curve  3,  After  Trim  Tank  Level  vs.  Time 
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Figure   IV-13 
Fairwater   Planes  Angle   vs.    Time 
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Figure  IV- 14 
Stern  Planes  Angle  vs.  Time 
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Figure  IV-15 
Pitch  vs.  Time 
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Figure   IV-16 
Depth   vs.    Time 
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v-   DISCUS SIONx  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   .  DISCUSSION 

In  the  design  of  optimal  controllers  a  major  difficulty 
is  in  determining  a  weighting  matrix  that  produces  the 
desired  response.  With  the  use  of  computers  many  simulations 
can  be  made  and  the  weights  adjusted  to  obtain  the  response 
needed.  In  this  thesis  no  attempt  was  made  to  weight  the 
velocities  with  the  errors.  This  is  an  option  that  could 
give  more  control  over  tne  output  response,  but  this  would 
also  greatly  increase  the  difficulty  in  selecting  the 
weighting  matrices. 

The  optimal  controller  designed  gave  excellent  results 
and  produced  steady  state  errors  that  were  smaller  than  can 
be  normally  measured.  The  magnitude  of  the  error  achieved 
is  a  function  of  the  weighting  matrix  selected  and  the 
forcing  inputs  to  the  system. 

The  depth  controller  was  not  designed  as  a  depth 
changer,  but  as  a  regulator.  However  it  did  function  for 
depth  changing  inputs.  The  depth  changing  rate  was 
determined  by  a  saturation  limit  in  the  depth  and  pitch 
error  inputs  .  The  limiter  is  required  because  the  planes 
have  a  finite  capability  and  the  use  of  optimal  control 
assumes  an  infinite  control  input  is  available.  The 
saturation  limit  can  be  adjusted  for  different  plane  angle 
responses  to  large  depth  changes.  This  would  change  the 
rate  of  depth  change  and  tne  pitch  angle  assumed  during  the 
transient.  The  magnitude  of  the  saturation  limit  chosen  was 
based  on  limiting  the  planes  angles  when  a  new  ordered  depth 
was  given  and  further  study  of  any  range  of  limits  was  not 
made. 
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Trim  control  was  accomplished  using  a  model  comparison 
technique.  The  results  of  this  technique  were  in  some 
respects  better  than  expected.  The  trim  controller  was 
accurate  in  the  final  results  and  its  corrective  action  was 
achieved  earlier  than  a  manual  operator  could  accomplish 
because  the  error  was  detected  sooner.  The  model  was  run  in 
two  degrees  of  freedom.  If  the  model  is  left  at  two  degrees 
of  freedom  something  must  be  done  to  compensate  for  the 
effects  observed  in  a  roll.  The  system  would  not  reliably 
detect  a  change  in  trim  in  a  turn  and  its  ability  to  detect 
a  change  in  trim  during  a  large  vertical  transient  is 
questionable.  However,  the  controller  would  not  mistake  a 
vertical  transient  as  a  demand  or  need  for  ballast  change. 
Because  of  this  last  result  it  was  not  necessary  to 
compensate  or  change  the  system  in  a  dynamic  vertical 
situation.  The  trim  controller  was  not  based  on  the  depth 
controller  used  so  it  would  conceivably  work  with  any  depth 
controller  as  long  as  the  planes  were  used  to  compensate  for 
pitch  and  depth  errors.  Filters  would  be  required  to  give 
continuous  average  plane  angles. 

Using  the  trim  controller  to  adjust  the  submarine  to  a 
"heavy"  or  "light"  condition  was  not  part  of  its  intended 
function  but  the  controller  was  able  to  accomplish  this  with 
reasonable  accuracy.  The  accuracy  was  largely  determined  by 
the  dead  zone  in  the  signal  channel.  The  deadzone  would  have 
to  be  large  enough  not  to  mistake  noise  or  ordered  depth 
transients  for  a  trim  error.  This  contraint  imposes  a  limit 
on  the  accuracy  that  could  be  achieved.  The  priority 
assigned  in  adjusting  trim  was  arbitrarily  assigned  to  first 
correct  an  overall  weight  problem  then  adjust  for  a  fore  and 
aft  trim. 
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B.    CONCLUSIONS 

Optima.1  control  techniques  can  be  used  to  design  a 
depth  controller  that  achieves  a  steady  state  condition  as 
accurately  as  desired.  Adjustment  of  the  weighting  matrices 
would  emphasize  those  errors  that  are  most  important  and 
also  adjust  the  amount  of  response  ny  the  planes  to 
accomplish  the  control.  The  insertion  of  nonlinearities  (in 
this  case  a  saturation  limit)  would  allow  the  controller  to 
give  a  reasonable  response  to  any  input  and  provide  for  a 
depth  changing  capability  as  well  as  the  depth  maintaining 
capability  designed  for. 

Steady  state  gains  were  used  instead  of  time  varying 
gains,  and  this  did  not  degrade  the  functioning  of  the 
controller  but  did  provide  a  significant  simplification  in 
the  controller.  A  controller  designed  with  the  methods 
described  would  provide  a  vast  improvement  over  manual 
operation.  The  gains  provided  are  only  a  function  of  speed 
which  makes  realization  of  the  controller  relatively  simple. 
The  controller  output  itself  is  also  a  function  of  pitch 
rate  and  depth  rate.  This  provides  an  anticipatory  feature 
to  the  controller. 

In  the  past  trim  control  has  been  purposely  left  to 
manual  control.  This  was  for  safety  reasons  and  because  the 
detection  of  an  out  of  trim  condition  was  a  matter  of 
judgement  based  on  the  observation  of  plane  positions  and 
therefore  only  good  in  the  steady  state.  The  passive 
override  mode  discussed  in  previous  sections  removes  the 
safety  constraints  in  the  use  of  an  automatic  trim 
controller.  The  simulations  run  demonstrate  that  detection 
of  trim  error  can  be  accomplished  accurately  in  a  dynamic 
situation.  The  method  still  uses  the  planes  positions  for 
its  determination  but  the  lifting  forces  of  both  sets  of 
planes  and  the  turning  moments  of  both  sets   of   planes   are 
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used  in  the  determination  of  trim.  Traditionally  only  the 
lifting  forces  of  the  fairwater  planes  and  the  turning 
moment  of  the  stern  planes  were  considered  for  the  trim 
error  estimate.  The  controller  made  accurate  estimates  as 
seen  in  fig.  IV-9  and  12  where  it  was  noted  that  the 
controller  never  had  to  go  back,  and  correct  an  adjustment 
once  it  was  made. 

The  additional  feature  in  the  trim  controller  of  being 
able  to  deliberately  order  an  out  of  trim  condition  provides 
for  more  positive  control  in  the  trim  of  the  submarine. 
This  control  is  a  function  of  the  actual  submarine 
conditions  and  not  of  the  assumed  conditions.  If  the 
conditions  change  during  the  submarine  operation  such  as  an 
ascent  the  trim  controller  would  automatically  maintain  the 
out  of  trim  condition  ordered.  A  difficult  feat  for  a 
manual  operator. 

C.    RECOMMENDATIONS  FOR  FURTHER  WORK 

There  are  six  specific  areas  that  need  further  study 
and   work.   They  are  enumerated  in  the  following  paragraphs. 

,  1.  The  selection  of  the  saturation  limits  in  the  error 
channel  was  arbitrary  and  the  only  criteria  was  to  prevent 
the  planes  angles  from  becoming  excessive  in  a  depth 
changing  situation.  As  in  the  actual  controller  there  must 
be  some  optimum  limit  that  would  meet  a  specific  criteria  in 
changing  depth.  this  limit  could  be  found  using  a  function 
minimization  program  where  numerous  simulation  runs  are  made 
and  compared  to  determine  a  cost  based  on  meeting  the 
criteria  specified.  A  minimum  cost  would  define  a  match  to 
the  criteria  specified. 

2.  To  achieve  a  zero  steady  state  error  and  at  the  same 
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time  make  the  controller  an  optimal  depth  changing 
controller  a  method  of  decoupling  is  needed.  The  methods  of 
steady  state  state  decoupling  through  the  use  of  cascade 
compensators  described  in  ref.  7  might  provide  such  a 
controller  with  a  function  minimization  procedure  to  meet  a 
specified  criteria.  This  approach  would  allow  for  zero 
steady  state  error  in  pitch  and  depth  under  most  conditions. 

3.  the  controllers  in  this  thesis  were  designed  for 
continuous  time.  To  provide  for  use  in  a  time  shared 
computer  environment  a  discrete  version  of  the  controller  is 
required.  For  the  trim  controller  this  would  mean  a 
discrete  model  used  for  comparison  with  the  submarine  and  a 
sampling  of  the  submarine  parameters. 

4.  The  control  inputs  considered  for  the  depth 
controller  were  the  fairwater  and  stern  planes.  The  control 
inputs  could  be  expanded  to  include  a  vertical  forcing 
function  and  a  rotational  forcing  function  around  the  Y 
axis.  These  inputs  would  of  course  be  the  shifting  of 
ballast  in  the  trim  tanks.  In  this  case  the  planes  angles 
would  be  included  as  states.  The  purpose  of  this  scheme  is 
to  control  the  planes  to  a  zero  angle  by  the  shifting  of 
Dallast. 

5.  The  trim  controller  does  not  work  in  a  turn.  A 
suggestion  to  correct  this  would  be  to  increase  tne  degrees 
of  freedom  of  the  model.  The  model  would  have  to  be 
developed  and  tested. 

6.  Finally,  the  options  allowed  for  the  pump  were  very 
limited  and  for  example  if  the  submarine  is  "heavy"  forward 
and  at  the  same  time  "heavy"  overall  the  designed  controller 
would  pump  ballast  from  auxiliary  tank  then  pump  ballast 
from  forward  to  aft.  A  more  efficient  maneuver  would  be  to 
pump  ballast  form  forward  trim  to  sea.   Even  for  only   three 
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variable   ballast   tanks  the  pumping  operation  would  become 
significantly  more  complicated. 
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EQUATIONS  OF  MOTION 

The  following  set  of  equations  are  referred  to  a  body 
fixed  system  of  axes  which  are  coincident  with  the  prin- 
cipal axes  of  inertia  of  the  body.   The  origin  of  this 
axis-system  is  located  at  the  assumed  center  of  mas3  of 
the  bodv 

Equation  of  Motion  Along  the  Body  Axis  System  x-Axis 

m{u  -  vr  +  wq)=  %  t*     [x     '  q2  +  X   '  r2  +  X     •  rpl 
n    Z  L   qq   n     rr         rp    rJ 

+  v  /3  f  X  .  '  u  +  X     '  vr  +  X     '  wql 
2    L   u         vr         wq    nJ 

+  £  i*   \x         '  u2  +  X    '  v2  +  X    •  w2l 
Z         L   uu         vv         ww     J 

+  -f  i2u2  [x  ,     '  6  2  +  X  .  ,   '  6  2  +  X  ....  '  6.  2] 
2      L   6r6r    r      6s6s    a      6b6b    b  J 


+  -f  Lz  X    ,  '  (n«  -  1)  v2 
Z  vvn' 


+  -5l2  X     ,  '  (n«  -  1)  w2 
Z  wwn 


+  %  X2u2  X      ,  •  (n»  -  1)  6  2 
2       6s5sn'  s 


+  -?iV  X      ,  '  (n1  -  1)  6  2 
Z  6r6rn'  r 


£W.  sin  8 


+  <Fx>P 
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Equation  of  Motion  Along  the  Body  Axis  System  y-Axi.s 
m(v  -  wp  +  ur)  =  |/*  [Y  '  '  f  +  Y  .  »  p  1 

+  f<4   [Ypq  'w+Yp|P|    'pIpI] 


+  T^     [Yv    '  *  +  Ywp    ,wP+Yv|r|^,Kv^w^|    |r|] 


+  ^Z3 


[Yr    'Ur+Y|r|6r    'u|r|6r+Y        'up] 


+  -?X3    Y       ,     •  (n«  -   1)  ur 

2  rn1  ' 


+  |i2    [y,    •  u'+  Yv    'uvt  Yy|v|     -vltv't  w!)J|] 


♦f'V*l,    '»r 


+  -£/2u2  Y  .      ,    •  (n«  -  1)  6 
2  6rn'  ] 


+  f  I2  Y       ,    '  (n«  -   1)  uv 
2  vn' 


+  2/2Yv|v|n'     '  (n'   "    l>v|(v»  +  w2)M 


+  fx2      Y  t  wv# 

£  wv 


+  -£x2(F)              v2  +  w2  (-w)     Sin  Wt 

2  y    vs 


+  £W  j  sin  (J)  cos  0 


§  Multiplied  by 

u 

U 

for  large  angles  of 
attack  near   -90° 
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Equation    of    Motion    Along    the    Body    Axis    System    z-Axts 


m(w  -  uq  +  vp)  s  |i*  Z   i    '  q 


+  £/*     [Zrr    T^  +  Zrp    -rpl] 

Note    1 

0,3     [z    .     '**Z  '  vri  +  Z  •  vp  +  AZ  '  vprf/1 

+  -y-X       L      w  vr  vp  r       u      vp  ^    J 


♦  fi3     [Zq    '«!♦*,,,,.    •«|q|«.*2w|q|'=|(v»*w»)»|H] 

+   |/3     Zqn«     ,(n'   "    »>«* 

+  |i2    [z,    'uz+Zw    '«w  +  Zw|w|     «w|(vJ  +  w2)*|] 

+  ix2  rzH,uiH +  zww"  iw(vz+ w^i  +  y»»»] 

♦  fiV    [Z6a     '6s  +  Z5b    '6b] 

+  f  X2    [z  ,    «  (n*  -  1)  uw  +  Z      .     |    ,    «  (n«  -  J)  w|{v2  +  w2)£|  | 

2  L       wn'  w|w|n'  '  'J 

+  |xV  Z   .      ,     «  (n«  -   1)  6 
2  6sn'  a 

+   -£x2(F    )           v*  +  wz  v   Sin  cot 

2  z     VS     g— 

+    £W.  cos  9  cos  $  a  Multiplied  by 

«  u 

for  large  angles  of 
attack  near  -90° 

Note   1 
when  not  multiplied  by 
u   add  to  Z      ' 
TJ  VP 
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Equation  of  Motion  About  the  Body  Axi3  System  x-Axls 

I^+(Iz-Iy)qr  =fz*  [K.  'p  +  Kqr  -qr  +  K.  'r  +  Kp,p|  lp,p|-j 

_ 

+  i1'  [K  p       Up  +  K  r    '  ur  +  K  v    '  *  +  K  wp    '  **] 

+  -£i3  [K,    '  u*+  Ky    '  uv  +  Kvjy|     J  v|(v*   +  w2)*l] 

+  -£  i3  K  •  vw 

2  vw 


♦!*"**«,    '6r 


+    Bz  „     ein  <{)    cos  9 


B 


10 


rr 


Equation    of    Motion    About    the    Body    Axis    System    y-Axls 

Note   1 
Iyq  ♦  0X  -/,)  rp    =  £  **     [M  *    '  q  +  M  '   r*#  +  M  fp    •  rp  +  AM  fp    •  rp*] 

♦   |i*     [M         'uq+Mlql63     'u|q|6sf  M   |w|q    «    |  fV  +    w^|q] 

^ 

+  •?!*     fM    .     '  w+  M  ■  vr#  +  M  •  vp#l 

2  L      w  vr  vp  J 

+   fX*Mqn'     ,(n'   -    ^ 

+  |i3     [m,     'u!+M  '  uw  +  M      |       ,      '  w|(va+war|'l 

2  L       *  w  w |  w I  '         •  'J 


+   -£x3     Tm,      /ulwl  fM        '     !w  (v2  +  w2)'|   +  M      'v2?l 
2  L       w|  \vw      I  '    •  vv  j 

+  f'3     "2    [M6s    '6s+M6b    '6b] 


+   ^  f   M        ,     '  (n1  -   1)  uw 
2  wn' 


+   iL*   Mw|w|n>     *(n'   -    Dw|(v2  +  W^| 


+  fx3  u2M  '  (n«  -   1)  6, 

/  6sn  i 


0  Multiply  by    rr  for 
large  angles  ofattack 


+    Bz        sinO  ncar  "90 

Note   1 

when  not  multiplied 

by    u  add  to  M      ■ 

-    £W  .  x..  cos  9cos  *  ,  U  rP 

i     ti  T 
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liquation  of  Motion  About  the  Body  Axis  System  z-Axls 


»-.*♦  Vx*)pq  =f/5  [Nf  ,f  +  Npq  'p^  +  Np  'p] 


♦  f  *4     [Nr    •  ur  +  Njrj6r    'u|r|«r+N|v|r«  |<v'  +  w2)*|r] 


*     [n        '  up  +  N  .     '  v  +  N  w       «  wp"l 


2  p 


+  -f  X*     N       ,     •  (n'  -   1)  ur 
2  rn' 


+  |/3     [N„    «u2+Nv    'uv+N.j     '  v|(v2  +  w2)*  |] 


+  $  l?     u2  N  .       "6 
2  6r  r 


+  -?X3     u2  N  .       ,     '  (n'  -    1)  6 
2  6rn'  j 


+  -?X3     N       ,     •  (n«  -    1)  uv 

2  vn1 


+  f'3  n„i„i„.  '  <n'  -  »)  vl(v*  +  w2)* 


v   v   n' 


+   -?X3        N  '       wv# 

2  wv 


+    £W.  x.  cos  9  sin  d) 
l     ti 


#  Multiply  by  —  for  large 
angles  of  attack  near   -90° 
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AUXILIARY    EQUATIONS 

<t>        =        p  ♦  i  sin9 

9         x        (q    .    i>cosQ  sinq)  )  /  cos<J> 

ty         =        (r+8sin<p)/cos9cos4) 

xo      =       uco88cos(f>  +  v(sin<f>sin9co30-cos4)sin(/j) 
+  w  (  8in<p   sin  #  +  cos<£>   sin8  cos||)) 

y0      =        ucos9sin0  +  v(cos(J)cos,/>  +  sin4>sin6siniM 
+  w  (  cos(f)  sin  8  sin  0  -   sin<$  cos  # ) 

z0      =        -  u  Bin  0  +  v  cos  0  sin({)  +  w  cos  9  cos  (() 

i 

U        =        (ua   +  va  +  w2)^ 

<Fx)p         =  fi'u^a^Vn'Wn"]  whcn  ^    <  r1 

=  ^  X2  u2  [  b£  +  V  n*  +  b3'  a'2  ]  when  k,   <   n'  <   kj 

=  ^  X2  u2  [  cx*  +  c^,'  n'  +  c3'  n'2  ]  when  k3  <   n'  <   kg 

=  I  X2  u2  [  dx'  +  cL,'  n'  +  d3'  n'2  ]  when  n'  <   lo, 


»    •  « 

i   t  1 

^1»  ^2>  k^       Sets  of  non-dl-.ensiot.al  coefficients  'ised  In  the  pro- 

,    f  pulsion  equation  above.   The  set  which  will  be  in  effect 

cl»  c?'  cl       a^  Qny  time  during  a  sirvila\.ei  r.aneiver  vlLL  depend  on  the 

t  ,       value  of  n  and  the  numbers  k.,V.p,k-. 

dl,  do,  d^ 
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NOMENCLATURE 


All  symbols  used  In  the  equations  of  motion  and  In  the 
auxiliary  equations  and  relationships  which  appear  In  this 
report  are  defined  below.   Any  dimensions  involved  will  be 
consistent  with  the  foot-pound-second  system  of  units.   All 
angles  are  in  degrees.   The  Fortran  variables  corresponding 
to  these  symbols  are  shown  In  Appendix  B 


SYMBOL  DEFINITION 

•  A  dot  over  any  symbol  signi- 

fies differentiation  with  res- 
pect to  time. 

B  Buoyancy  force  which  Is  posi- 

tive upwards . 

m  Mass  of  the  submarine  including 

the  water  in  the  free  Hio^rig 
spaces . 

j^  Overall  length  of  the  submarine 

U  Linear  velocity  of  origin  of 

body  axes  relative  to  an  earth- 
fixed  axis  system. 

u  Component  of  U  along  the  body 

x-axis . 

V  Component  of  U  along  the  body 

y-axis . 

« 

v  Component    of    U    along    the    body 

z-axls . 


/ 
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Command  speed:   A  steady  value 
of  u  for  a  given  propeller  rpm 
whenc^.p  and  control  surface 
angles  are  zero.   Sign  changes 
with  propeller  reversal. 

Longitudinal  axis  of  the  body 
fixed  coordinate  axis  system. 

Transverse  axis  of  the  body 
fixed  coordinate  axis  system. 

Vertical  axis  of  the  body  fixed   ^ 
coordinate  axis  system. 


Distance  along  the  xQ  axis  of  an 
earth-fixed  axis  system. 


Distance  along  the  y   axis  of  an 

7  o 

earth-fixed  axis  system. 


Distance  along  the  z   axis  of  an 
earth-fixed  axis  system. 


Component  of  angular  velocity 
about  the  body  fixed  x-axis. 


Component  of  angular  velocity 
about  the  body  fixed  y-axis. 


Component  of  angular  velocity 
about  the  body  fixed  z-axis. 


The  z  coordinate  of  the  center 
of  buoyance  (CB)  of  the  subma- 
rine. 
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Angle  of  attack. 


p 

Sb, 

Dt 

Sr, 

Or 

Ss  j 

os 

Anfclc  of  drift. 

Deflection  of  bowplane  (  or 
•  al lp  lane  ) 

Deflection  of  rudder. 

Deflection  of  sternplane. 


The  ratio  u  /u. 

c 


/ 


e 

c 


Pitch  angle. 

Yaw  angle. 

Roll  angle . 

Mass  density  of  sea  water. 

Weight  of  water  blown  from  a 
particular  ballast  tank  ident- 
ified by  the  integer  assigned 
to  the  index  1. 


v 


V 


0> 


Angular    velocity 


V 


ti 


Time  . 

Location  along  the  body  x-axts 
of  the  center  of  mass  of  th$  itn 

ballast  tank  when  this  tank  is 
filled  with  sea  water. 


Ill 


(Vp 


Propulsion  force  (see  auxi- 
liary equations  and  relation- 
ships) . 


y 


Monient  of  Inertia  of  a  sub- 
marine about  the  /-axis. 

Moment  of  Inertia  of  a  sub- 
marine about  the  y-axls. 

Moment  of  Inertia  of  a  sub- 
marine about  the  z-axls. 


P'    KP'     KP|P|    'Kqr- 
r*    Kr*    Kv    kwp«    k*. 
KV    Kv  jv|  »    Kvw    KSr' 


Non-dlmens lona 1  constants  each 
of  which  Is  assigned  to  a  parti- 
cular force  term  In  the  equation 
of  motion  about  the  body  x-axls. 


M.  '  ,  M   •  ,  M   '  ,  AM   '  ,  M  ',  M,   ,   ' 

q  *   rr  '   rp       rp     q    |  q  |  6s 

u.  .    M  '   M   '   M   '   M   ',  M.',    Non-dlmens  lona  1  constants  each 
I  w  J  q  *   v  *   vr  *   vp     qn     * 

of  which  is  assigned  to  a  partl- 

M  '   \{     ,      Mi   i  '   M   '.  M   '.  M ,  ',  cular  force  term  In  the  equation 
w  *   w|wj*    |w|   '   ww  *   vv     6s 

of  motion  about  the  body  y-axis. 

M«b''  Mwn«'»  Mw|w|n'''  M6sn'' 
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r    pq     P     r  '   |  r  J  6  r     |v|r 


p  '   v     wp     rn      *     v 


Non-d inens lona I  constants  each 
of  which  Is  assigned  to  a  parti- 
cular force  tern  In  the  equation 
Nv|v|  '  N6r  '  N5rn'  '  Nvn*  '  Nv|v!n'  *  of  motion  about  the  body  i-axls. 


N   ' 
wv 


A      •   A      ■   A      •   A  •    t   A      •   A      « 

qq     rr     rp     u     vr     wq  ' 

X   •   X   '   X   '   X     '   X     '    Non-dimensional  constants  each 
uu     vv     ww  '   6  t6  r  '   6  s6  s  * 

of  which  is  assigned  to  a  pnrtl- 

X,  ,  *   X    '.X    ,'   X       '     cular  force  term  In  the  equation 
0  b6  b     vvn     wvn   '   <5s<5sn' 

of  motion  along  the  body  x-axls. 

^rirn'' 


V«  Yp'»  %<,'«  YP|p|'  V*  Ywp'» 


Yv|r|''  V'  Y|r|*r''  V'  Yrn''» 
V-  V'  Yv|vf'  Y6r''  Y6rn'' 


Yvn'''  Yv|v|n,,«  Ywv'  »  (Vv» 


,   Non-dimensional  constants  each 
of  which  is  assigned  to  a  parti- 
cular force  term  in  the  equatlo 
of  motion  along  the  body  y-axls 


V'  Zrr'«  Zrp''  V  '  Zvr»  Zvp'' 

AV'  Zq''  Z|q|6s'  Zw|q|'' 
Zqn"  **'•  V'  Zw|w|'«  Zjw|'' 

Zvw' '  Zvv' '  Z6s' '  Z6b' »  Zwn'  ' 


vvn     6  s  n      z 


Non-dimensional  constants  each 
of  which  is  assigned  to  a  parti- 
cular forci>  term  in  the  equatlo- 
of  motion  along  the  body  z-axls 


113 


oo 

r«- 

f\j 

0 

0 

in    0 

in 

a 

0    m 

ro     0 

»-< 

UJ 

•    r- 

ino       • 

0 

oo  uj 

00 

1     0 

voo      1  in 

• 

QdZ      Q 

aim 

11  mo 

0  •      11  i-n 
•  1      i»o 

11 

UJ>-4         O 

UJO 

a:(M  • 

> 

Y-CC       OO 

>--. 

oo  1 

ti  11      ico 

> 

LLXOQ 

UJO 

cto  11 

>CD         »    • 

2I<M 

SSI  »- 

CJ>       CCQ 

QO-^M 

>o      Is-  II 

»o 

m 

<CO       O 

Z             • 

x  .0:0 

^J^     o± 

ino 

0 

a:3az 

KH          O0   II 

c  II    —lO 

•-  ►     1  > 
coin     my! 

Oo 

0 

<O0UJ>-< 

Z         "-Hif) 

(\JCL>0 

0   • 

• 

Q.      H->-l- 

0        <t 

o.r>-   • 

00     0  »• 

•  l 

1 

OoOLLoi 

•-"UJUJ    •> 

OX   ~ll 

•0       »in 

l   11 

11 

LUZ>— *•— •  _J 

i—X-Jin 

•   ►roO 

1      •       OOfM 

II  h- 

oc 

IMJOQ. 

—*\—  Oo 

11  ooa. 

M    1           IO 

irp 

0 

H-t-     uu 

v>      Zo 

300> 

— 1  II          II  O 

~«o 

z 

<xc*:cli— 

oz<o 

20    •   - 

:<oo     o-o 

3C2 

► 

<^ 

^ 

f-     •  _JtUl/}Z 

O-UJ       0 

XO  II  <\J 

f— 0     <— »  • 

^2: 

m 

—J 

- 

OSDQ      < 

TUU    • 

►00:0 

12>sl      Cl  1 

2    - 

0 

oosicxrx- 

>2Z  I 

in  «>o 

••in     i<r- 

►in 

0 

IQmOcidw 

CD       <C  II 

v0  1    »o 

mm     >o 

0 

ooujoo         Z 

_K\In}-0  II  s0    • 

>— 1-*       »o 

00     0 

• 

QiUJ       LUQO 

O'/kkooqio  it 

—<OOvOO 

00      1 

11 

ujarc^xujo 

lulu        ••  1 

•  rf    »CC 

oooo> 

♦  00  uu 

y- 

>u.ni--j 

XUJZ  — 4UJ 

1  II  X  1  O 

•   »o  1  ^i 

ii    .00 

0 

O       U_       _»LU 

ooc£a:or^>  ►  n  > 

1      1    OUJ    »• 

:*  i  0  • 

0 

u_     zox 

►-•CJMJJO-O>(\j>    - 

11  11    «om 

2;  II  ON 

> 

uuoaj—<cc:K 

_jujh-   »ooxo-irf| 

'JioC  1     «o 

-±   *  1 

Z 

Z      O      1- 

O-Qoo  |   vo 

1  *-o>o 

r-sjCX  11  r—  | 

incc  11  11 

•* 

<^MtUO> 

2:            11     « 

^0-!0 

»rf>JCL  1  UJ 

02CXH- 

O 

_JLLO<X)OCiD    « 

'Oouj»-*mN  «>o 

-H     "ijl    II   O 

1    »->o 

O 

O.UJO                 QONK  II  O  |     ►    • 

OCOrvlH-    • 

ujvj-2:o 

1 

Cti       hLUDZ 

:o  -s  1—  *.t«JO  11  — ■  1 

00  -o-t 

00  -a. 

UJ 

in 

_JCJ>O0O0q;Uj3<             lA"-! 

•  O™  II 

OOOO  1 

•  ooz 

O 

O 

• 

<ujo:z5<i:2c 

1      QoOsf    * 

■  \   OOh- 

•    •    *CC  II 

xj-        •        •       •• 

• 

• 

11 

oQiur     t—UjooH-ujvrvt  11  x  »a 

1    1  O^ac 

11   1  o<t 

r\J 

0 

00 

►-h     on     wxi/i 

!►-«     UJr-00  ••  1  0 

11  11  11   »>k 

00   II     II   O 

l 

11 

_j 

j-voDU-uja: 

Oaico  I 

-so  11  a. 

000^3—*    •> 

2(*-H    I 

II 

*- 

IU 

ho:         u.Ziua:(-t-ooujxM>> 

fvi  >  j£  cd  in 

►>3:uj 

a 

O 

0 

*ujzo<<t-i.uo     Luor— 

«o>  •> 

»rvlrvlOO  ^tinSI-HO 

z 

O 

► 

>>—0  1   _JUjCj-0>-0    »>o 

■— *o  •■o 

f—    •"  *-0  1  O— *    "S!    • 

•» 

1-^ 

T-4 

2:        a:<o-o 

X_J        II  CO 

—4     »LPlO 

OJ3  0    »JJ  1 

OCM    --J- 

in 

no 

0 

<zooo<r        i>i^joj-ooino 

00   »  1  in ujoo m  1  <j-  >^ 

*.    • 

• 

o£Oujz<ujooz 

£*•<    »2    « 

•  II  00 

•  0011    •moooii  — 

lOM 

00 

II 

Omk         x«-«<uji— tin  ••ion  cQOO 

1      ♦  il  CCA    < 

>  •  «ooooo 

•  II 

1— 

Oh-  <UJO0h-       —1 

>>i~  1  m  11 

mo  •   • 

II    1  — <0  |   (NJ  ll    1   OQ.    « 

•  «o 

ooc: 

-J 

C£<_JXUJ       1— QO^-<       rO>>cDII    II 

a.  11  2^  11  11 

a.  it    »z  11 

1  om 

11  0 

UJ 

a._io>-_i  «z<r 

P"oo   •»— 1 

XOQ-h- 

>~*r-*  v -4.Q.a: m  n    •-> 

■  11    »o 

»—  ► 

0 

—j  2!     co  uj  lu  ai 

UJ«— iUJlPi    » 

■  ►y^n 

^^JC3^^sJsO  it  S 

■  ■ST.— icd-^-^a:  l    1 

no 

*» 

S-S«J32 

Z        XfJ    OLTi    «->-0 

-— •   ►O   ►¥ 

:  ^Qo^ 

•—I  II  UJ 

0   • 

0 

(-4  00  0.-._JLLI.^<[00<.J  11  0 

,_<LA     «-Qi 

m^c\!  1  o>   »--j---<2  1  rH 

i>2.rr> 

00 

• 

/K 

OO        ttccO'"1 

_IUJ<0    1 

oro--> 

NnJOUJO'AOS    ►luZh2-< 

ct  it 

0 

O 

*         oou_oo 

ClZLUDUJOOh    « 

0  -OOOOO  -ino  -ZZ^rvi  -co 

in 

O 

tu  *-*     >     2: 00 

'     <a:  •"J-oooMT.omo  «o  i 

oincNj  fA 

1    ••   »^i 

\OQ 

00 

r\JXX 

0     z     q     2:     — 

-J_j     0^ 

•  t  »o<: 

1    »^-    «rn    «ll 

l   •r\ioinr^mu"\>— •  i 

:\j  •  •» 

ll 

■>^~ 

II— 1  »■ 

hX 

Q      hhUUUJLUOU- 

<ra_uj  •— • 

1   1   1  oc 

>  1  O  1    1    1  c 

1  1  00  1  00— <o. 

—iOO     —^2 

a  1 

1 

1   1  00 

CCZZIwOhU       -Jinr-0  || 

'  II   O   !l     II     II       « 

no  •  11  000  • 

•  II      »V~  II 

•— < 

L)N>-f 

-J-J  •  • 

_JH-       <I "— >  y-*  h-  H-       0l 

•— < 2TO  — *    • 

1— 00H-    *   1 

t—    •!— >— >fs-l—    •   1   h-    « 

1     •  *H 

II  hoOl-l- 

-J1— 

l^-(  h 

-200 

003      2Zu:<I      OUJ—J 

1—  a^z^-f-oao  I  II 

O  1   03-*  II  O  1    II  O  1 

1    II    II  h-O  II  JOZ 

II  II 

II    ! 

II  II  II  II 

00-        CQ<|--     •2Th-Q.Q.UJ<ll    II 

o^r_>  u  >o  11  Go>  Jio  11  000  11 

II   >QJQiOOOQ.O»-« 

fNJXX; 

>-<cin_> 

Z      3S23ZO< 

uj|—     oxDo>xssaoa-«<Qa3ua:>ii:2Nj<QZQ.a. 

UN>N<« 

o>-<    wcQOOu.a;iiJ 

>ooX-J>-« 

xx>>>-^r^i*v^^^: 

:zzzzzzzaa.>- 

:     -j 

'-J  — 

1 1— 1 4- 1  <x  <l  <t 

UJ    •            3      hh            :£ 

0 

C£ 

_j< 

X  -J       LU  00  UJ  I—  UJ  UJ  UJ  UJ  LU 1—  5"  21 

2: 2: 2"  2:2 

2:3:2:22:2 

22:222:2 

2ZZUJI— Ct>- 

1 

UJOO      _J      (XO^IlTI,-,<<''<<<<«'t<t<'<<I<<<<<^<<<OOOi/)(-(- 

Q      H- <.<  2  3 H- 31  h- H-  Q, a.  u£ 

acccu^a 

c^XcccLCX-alaLaiCCcccc 

aa^aouHZZw 

N,^      <-> 

•:r<i<t<x<i<t<r<x<i<r<j<t< 

<t<i<<r<<r<'--K^:zzooz 

V^V.-K-h-*!-  tt  ■«-  tt  #  tt 

■>>  #  ^'Xana-aaa 

Ci.Q.a.o-^i.Q.cxQ.cxcLa.a.ci.Q.CL'— i.-«»-ioo— < 

cD 


11^ 


_i  _« _<  ,-«  — i  .-i  r\j  cm  r\j  r\j  r»  i  r\i  m  m 

o  -^  <m  ro  >j-  in  <o  —4  cm  m  vj-  it*  vo  — *  .m 

o<K"<r<r<r<i<<i<<r<<<r< 

••<  <i<:  <r  <<  <i<i<  <r  <r  <  <  < 

<j    <^>  <~>        <^>  o  o  <<3<<<<i<<<<<r<<i<(< 

k—         _)         _J  »—         _J  _l(-  _J        h"  _J  h-<Tw—  ^w—  —  «-www-.-—^w 

O       v.       V.  O       *>»  -5fO  *CD  #  0<LULXJUJiJUUJUJUJUJUJUJUJUJUJ|JJ 

O      H-      V-  OH-  h-Oh.  f-Q  ►-(-      Q-D^DDDDDDDDDD^b 

>  O      O  SO  OQ.O  OO  OO      oc:uj_j_j_)_j_j_j_j_j_j_j_j  j_j_j 

>  O      O  N      O  O      *£      O  O      ST  QQ      ZD<«««««J<« 

oooo  i  o>o>oo  »  o^ooo-o  t  oclooqto  i  ooqcod^o  i  _j>>>>>>>>>>>>>> 
•  •  •  •_)  •ii  »z  •  •— )  'Si  •  t>  »x  «^  •  »-vi  •  >-  •  •>-  »\i  »>j<i  ii  it  n  ii  it  ii  ii  ii  ii  ii  ii  ii  ii  ii 

OOOOSIO  I  O  I  002:0  I  OO  I  O— <0  »  OO  I  0>—  OO  I  O  I  OM><aoQUJUL<'I]UQiiJ!J.<CD 
II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    II    11    II    II    II    II    II    II    II    II    II    11  <<c<i<<i<icococQrorocaoo 

Quju.<ccoOiijj.<cDuaLUu.<iriOQUJU.<cDuQuju.<cooQuJu lu.a.u.u.u.u.u-u.u.u.LLu.L.-a. 

«<iQCQCDC0.t)CDjOUUUOQJQQQClUUJUJLUUJUJU.U.LL'XU.U.UJOGOG3OQOaOOQaO 
<<<<<<<!<<  <C<<<t<I<<  <<<<<.<<<<<<<<<<<<  OOOOOO'OOOOOOOOO 


115 


-^  £,    o  - 


a.    . 


, 


0    - 


•— 

<\jf^ 

— K- 

f\J(\JCO 

H-O  — 

eg 

OQI- 

OQ 

O'-'C 

1—  t-O 

OCLQ 

OOH 

c£  «-<ra: 

ooo 

►o>o 

(•5 

o 

r<~>  ro  ro  m  <t"  vt  >*•  vt  ■J-vt^iniAinininvOO-O'O^OO         OH 

<<i<i<<<<<<r <i<<r<K<r<<:<<<<K  •  ►ooooao  •  h—     — — 

<<<<<<i<r<<<r<<r<<t<r<<;<r<<r<r<  omooooo  »oo<     oo— » 

<<<<<<<<<!<<<<!<; <!<<:<<<<r<  ^rO>2Q.3o:Q-'   »c£        •   «r-^ 

«»*%^—  ~>^<-«^« ^-— -^» w,~^»_>,,_ -_. «  <»* 'w— '«_• •_•>«•  _j    #•.».».•.  ».«w. lOai^OO    • 

LULUUjmUJUJLULUUJUJUJLiJUJLUiJJ'-UUjUJUJUJUJUJ  QL ITV    •     •     •     •     »_JC£-*2;Q     •     »0 

Z5r)3or>r>D3DZ)r)r5r333-DrDr5Z)Z)Z)^       onooooo^ojuj  i  an-uj 

-J_J_J~J_J_J_J_l_J_J_J_J_l-J.J-J_i—l_J_J  — J -J  V--— ——Oh  c£OQLU_j    • 

<<<<<<<<<<<<  <r<t<<t<<<<t<<  ZJ-IJJ-IJI-ZO      a    •    «Q_ 

>>>>>>>>>>>>>>>>>>>>>>  -'alcccccc^mz^i-uj'^cjicc-sio^ 

II   ll  II  II   II   ii  ll   II  II  II  ll   li   il  II   li   II  if   il   ii  il   II  ll  II  (JOOOOO>-i  il  Z-JQujuji—  II   it  oo  —  I—  Or\i>—  H 

OOUUU_<CDOQUJU.<cQOQLULU<C0OQUJU.  Ill— I— I— I— t— I—  II  I'-'O  II  LOooi/121  a.  Q>— wOO      t/)a 

OOUOQQQQCQUJUjjjUJiUiuu.U.lIiiU.u.'II  hZZZ?ZZ_IOII  2a:OQQ^.5.  II    II    llOOQS 

LL  LL  LL  LL  UL  U_  LL  LL  U_  LL  LL  LL  LL  LL  U.  U.  Li.  U_  LL  U_  UL  LL  >  3l  —<  »-•  —•  ►""  — <  *-*  — 1 1—  2  <  UJ  —  —>  ~-  h~  N-  (_)  I—  <  III—  ^.|— 

OOOOOOUOOOODOOOOOOOOOOh      uj  II   il   II   II   ll   II  a—i<      (y)u_LLLLi/)oooOLOooi/)au_<y> 

<l-  z 

>ct  <  co         r- 

or</>  a. 

iuO 

QZ  -H- 


o 

eg 

1— 

a 

a: 

h- 

_j 

o 

a. 

to 

* 

— » 

—» 

i— 

— < 

LO 

• 

O 

3 

•    < 

LU 

ILICO 

• 

LUSIQ 

s: 

"z—  +  ~*-eo 

116 


> 

2 


< 


-J 

X 

-J 

u 

o 

»— 

a: 

1—4 

■«■ 

CL 

/"" 

z 

«•« 

•-H 

OO 

00 

O 

■K- 

o 

^ 

^ 

X 

^ 

o 

«-» 

1- 

—_) 

•— < 

_I-J 

CL 

-JO 

«^ 

Oc£ 

oo 

a:— 

o> 

O 

_-z-^ 

in 

o-< 

,z»-«x 

f— i 

#  -J 

M(/)0 

<S-<J)^JU 

>_JOO#-  V- 

^— • I  — I » 

+  0-K-  or«-H 

•V^V^^v 

-^OiOi  +  O. 

roarooco 

X-- 

(  -~w 

^QQO 

c\j(\jcsic\j 

O00 

-JZ 

X*  X  -K-  iJ>i~i<J>UU^_)^) 

l-O 

-J-J-- ' 

L>rr'^ocn_j_j_j_i_j_i_j 

IOwOJO^vO^OI-OQQSSSNwSS^ 

— o 

—o^ 

Qu-K-^O 

-«•  :scl-i 

Qi-JJ-  —  •>* 

-*  J2-» 

cus:  >-•<%] 

a:  I  oo# 

XOr  |  * 

+  &LJI-) 

-*— {\j3-.jr3 

OOft  -M-  ^tt 

•    ••H-0-^"0 

•-i— «n£3:rvl< 

«  ••■«■  X-M-  •*■ 

ooar  +  -h-o 

•   «a:a£3rr> 

OOX-*  ->  ■«■ 

•-  »  +  >2D 

OONft  "-£■«• 

•   •  •»■  ct  !x  r-j 

00 

—«—«•{«•  >  *  < 

2T 

— .  1     1   3Xr\J  + 

O 

3  •••-*+•>}■  ro 

►— i 

AC 

:o>aa«  -::- 

1— 

«■ 

>w*.Oft  >•!!• 

rn 

m 

<~ 

-» ?>-;2X>*  r) 

O 

UJ 

-j>sa.oa:  + 

•~«  oo  oo  —  -*  >  -fc- 

a. 

a. 

UJ>— 

«- ^*>t-ZZ4  _J>— « 

+ 

•v 

OCOOOOOOOO  OVIS- 

ocoiju5:>« 

Cu 
* 

— «  Za#  OQC-m-  CTv  j>>-<->  ->  -K-  alo-ocOGioOcoaLm 

lu  o— — od^t—  iMMarso^QOQQQnQ< 

31  mzIwwO    •  •-"-  *  -;:-  -a-  -::-  *-  ->  -a-  -:-  -»■  tt  O- 

•-h  l-^OooOnr-P^    03DDDDDDDD  + 

»-<  < OOI—  CD O <  WMT» O  -If  *•  tt  #■«■->  -if  -II-  #  -H-C\J 

»00  Z>#«-<0 -«•>-■>;■  -K- co c<  oo  cQ r>  O Z3 Ij  r5 D  r>< 

uji— q  iuQzjo.  -ft  a:  +  oocD— iQQQ*  ■«■  -tf-  •;<•  -:$■  ■**•  -m-  o. 

~«      5IoO  +  — iOo^lu  I  «3wQ.OQ  ||  n^oo'Qc^oocOciCoTcoa;  + 

II  00>-<OOf\l  II    112       II  oo  II    II   II    II    II  <QQOQOOQaQO-<-HiM-((MH      CO  rO  CO  CO  nO>OU_U l~"~»— 

SQI-wOO  2:0.— l>t-Ot~h-l-«'XXXX>NN^SXZ<CD(JQliJa.a:<<<<<  II  00  !|  ||  +  +  4.  + 
i.  II  II  HQOiSI-CCOOGOO^K'JII  II  II  II  II  II  'I  II  II  II  CL  0.0.0.0.0.  <  II  II  II  II  l|  £.  II  O.O.  II  II  II  II 
I—  OK-<  II  I—  ^-f—  Z<Q  It-  QQQOOh  —•  c\)  m -<  r\j '.ft  ~*  r\j  rr» .-«  II  n  11  11  n  11  uj>3Ql  G*oi  2  "S  2> -Hi'M^^- 
tOi^oo^oOOoOoo  Q^Q3  ^  <0  O0t00«<<  <aDOOOLUUJU.<cOOQ  uJU.ZaDcQ.Ocom  >>  X)cD<<  << 
QQQQQOQQOJN  +  Q.  >OC  QQQ-Cl  O,  Cl Cl  Q.O.O.  O.Q.  a 0.0.0.0-  Q_  Q- ►-><<<<<.  ><<<o000o0i/} 

-I 

z 
o 


a> 

<\J 

— <x 

r\) 

<M 

< 

117 


UJUJUJUJUJUJ 

QQQDQO 

U.LLU.U..U.UL. 

O  O  ft-  ft  ft  ft  ft  ft  . 

-J  _J  U_U_U_U-U_U_ 

""•^  "^  •  <XC0OQUJLL 

>«-«  --»  •  LL  LL  UL  U_  LL  Q. 

•  #">  oO  •  OOOOOO 

•  q.ivj~.  s;— >  t  oooooo 

•  -ft  ft-  D  O-ftX  ♦  4.4.4.4.4.4. 

~*  — <0<  JNU  •  UJUJUJUJUJUJ 

X  21  ft  4  "^-ftl—  -ft  rsOivjtv)rNjM->J 

O  I    Oh-  —  -»■•-«  —  #  #  *  *  •>>  # 

—  t—  0  4LL  — »OCL  t—  UJUJUJUJUJUJ 

0£»-<  #3+  X4-w  <  <CQOOUJu_ 

■r- a.  occi-  03:00  ft  u_a.u-u.u,u_ 

3«—  -ft  <■<.  h->o  in  oooooo 

#00  _j-*«—  «-<<o  •  oooooo 

-JO  sr—t-K-  cl  ■»{■•«■  0  444444 

SO  +3—  O             ~-0  —  ~t  oooooo 

I   -ft  CL— «_J  _J  Z-K-  h-  CM  rvJivslrvlrsjiNj^j 

CL  — »  3rvl_l  "v  i~«-H<r  J  ft -ft -ft  -ft  ft  ft 

-ft  _l  rQ-ji-O  — .(M  003-ft  *-  OOOOOO 

3_j  <oa:  — o  ft-— im  u_  <tccooujll 

ft  O  *  +  «~»       X  -J  CD-SI.  -ft  U.  LL  LL  U_  LL  LL 

-JQi  33:00     O^  N  +  C^  'A  oooooo 

S^-»  >>Oa:h--»  CQ3— 1  •  oooooo 

4-z:  <<  o  ft  ►-<  c\j  1  ^  fNi  m  4.44  +  4.4. 

Q-1-"  -ft  ft -ft  o a. -k-  021  r>-  oooooo 

>GO  03— »ft  —  -ft  O       -ft  >l—  — 1  tvJNjrNi^rNJM 

c04  co-ft  :r>-oooft      3<ll  *-»  ft  ft  ft  ft -ft  #- 

<o  <IHONO#-     >ft  ft  o      +  oooooo 

ft-  .J  ft  3K-'-<OgO(\]      <3in  _J      >  <CDOOUJu_ 

2V  Hr^M|    «^(J        ^H     .  "Vf\)-ft  LLLLLLLLLLLL 

>-*            c*3cl-^-^4._j     o-tm        -.00                                                   oooooo  -> 

<>     o     ■- <Nj^-a__J3v.a:.-H,-ir-         c£_Jftf\j  oooooo ■< 

ft -ft  -ft       2  +  oO*_Ht>*2SH           ft*v.>0  4444.4-4.  1K.X 

GCS.  ~.     rvi<\JO0-O>ft  — •-!+-           3>Z_J  CQCDCOCDCOCQ  I— CfO 

CQ>  U-*  +  4!-  uaD^-*>a.JT!M  I             >ftwvv  fvdrvI'vi^fvj'^J  Q.OH- 

<<  _JQ-r<:ft  4<I—  ^-^-H-  4-  ft  3:           <;.— 4  — .  -ft  ft  ft  ft  ft -ft  LUOO>-« 

ft  -ft  V.-ft  ft>Oft."?^>4.X>ft-i!-            -ft  3  O  -»  CO  CD  COcO  GO  rO  QQQ- 

.-<>  >Q->ft  OQ-i-H^otNJft  "?:0           a>— I— I                                                    -  <CO(_)OUJu_  •■   ••  •» 

CC-*  -ft  Qi-JS- >"^^00  4  _l»-i  — -ft  -ft             — (<-V_J<f<J-vt^--J-<t-  LLLLLLLLLLLL  UJUJUJ 

-<>  >rvja:  >  3  CV  ft  3^.40.33  0>ft  DO<OUOUJU.  OOOOOO  STSTS 

">—*  >-4">r^)->f-  -HCQ>-»Qf:-i!-  35*  ft  —«>  ■«•  n£  GO0O0O'  OoOoO  OOOOOO  mmh 

0>>-      4-Qirvj  4  o*:->j<ra:->  a.  "C^-rMQ.  "="--<— -»  +  •+•  ■+•  4-  4  4-  Z4-  +  44.44  f-»—  h- 

•ft  4  4.       a.ft+CM-£+X>-frftc£>440ft+>-£zrorOrorn<-cirO  0<<<<f<<  •.».»■ 

Cl<X>      -ft  c£CLft  3ac"  I  >c£QCSr\JO-J— -a-.-H->r  ^<cCOOujll  i-Hi^rvirsji-jrvj-vj  r-t-nro 

#■}<■■«■      Q.a:-»<-  «■  tsj-ft  a.  ->,:-  ^5"  ■*■  -ft  jvy*  2*0:00000000000000  ^- ft  ft  ft -«•  ft -c-  *•  *•  - 

023^>iM>T:4.0tt>+ J-^-i!-N~>i:  +  Zv  +  +  -!'  +  +  +  0<<l<<<<  r^fMro 

O-ft  ft  O  +—  ft  -ft  Oft  3-^  a_Q.-ft  >  O— •-'-ft  3  4-  -»c\ir\jrvjr\j  <^f\J  S<cDOQUJ(i.  w.^»>- 

>0-><a:-M- 0-3ft  rx-ft  >-:}•■«•  c<ft  ft —1  j  a.ft  Q.x<cqoouju_                            u_u_u_u_u.u.  OOO 

■ft  33  +  ft  a>2DOl'-<aCL>>DJ'3j:>-j  oooooo i/i^i^<cooOUJU.LLOOOaOO  333 

o>->htf*NNii-^3yyas>*oaz2ah  +  +  ■+•  +  +  +aao.aaaOouuuau  a:  a:  a: 

Jn»s.u.>-u>— o>-^-««s«-5:oa:z-zzH"4-<^^-*-'+  +  +  +  +  +     —  ^-^->—w^  ooo 
+  +  +  +wj++(sj  +  «  +  ws*N.^5:«-  +  ^-a<'flOQii;u.<cDuQiiiii^ii  11  11  11  11  11 

II  ll  11  I—  11  11  11   11   11   11   11   il   ll  ii   11   11   11  00  11   11   11   11  »^ </i tr) go  1/ )<•)'./} go <yi oooooo  1/) s: i—  1—  l—  t—  I—  1-  r_i 1 

•-•  cm r<-> <>*-.- •rxjmvt— <<\jro<t-— «c\jm-ro~ irvj^isroo  11  u  11  h  h  u  m  11  11  11  11  11  jqodooq  _j_j~j 

COCOCO'— CDOOOOQOQQUJUJLUUJOU.U_U_u.  0<rQOQUJU_<cQOOUJU_H-.OQ0QQQ  <<<l_J 

oOoO^Oft  00  OO  00  GO  OO  oO  00  00  00  OO  OO  00  00  -ft  GOGOGOGOO^OOGOGOGOoOrNlrvlr^^sl^fsjv—  0>30-0cx:  ^KJ^X. 

<  LU  ~Z. 

O  I  fc— « 

O  l-Q.  S 

U.I  0:2:  CL 

0<  UJ 

•«■  0O0O  I- 


118 


o 
o 


a: 
o 
u. 


o 


f\j 


o 

-J 


< 


2:       II  C\J 

~                                                                                                                     <       C\J  -* 

—    o                                                             z    2:  ► 

2:          •>                                                                                                             ►— 1         ••  — < 

-    w           o                                                    a:    cm-»  «-» 

_       CO            — — *                                                                                LU       5IsO  — 

CO      —           2:0                                                                    2TLU      —iro  -J               r\l                         O 

«*      O            •+-                                                                2IQ      21— «  •— »            »          -,          vt 

lu       ••        #-»                                                          •.•      —uj  Zv:-»        o        "^          •> 

z>    vo        2:0                                              •-"►<»  z*:_i       r\i         ►       o 

<       <OOCQ                                                                             -»        ••—.—«  2:            -J— ^           'T)2r       ^ 

>          »  ►  11  — ,                         •••••••     o-  -^.^w  •»         11  <r^         —1  ••    — o     — * 

^,^_<^.o  •  •  •  •  •  »oooooo-«     ^-  *-iimo  •-<        i<;— »^     uj<-»„j— ■»  <  11      m 

2:    o  11  11  00  •000000  ti  11  n  11  11  11  11  uj  •- *•>..-<  •»— ioj  it         *:i/)_    r>fvi  u  -3 11  — »    — 

o     1-1^005:0*11  11  11  11  11  11 — »— »—.-*#«»— »o ,o— .2r>o-— o  o_>        ^aD<r     z^-j  •— »->     -z 

»-i       002:5:    «*UJ .—.«-»'-— .2:  ^I  2: 2:^^:2:  2 «- 1—       «-»-;Z  Q                        <w!^>-i   1         _j-j    »qz 

»—      Z          — 1   •  —<r\jcOvj-in  0 21 21 2: 2! 2: 2: 21  ►-' aj < o  1x1  < >-•  »-toO(\j  1  oo^kj^w  »mq  | 
O       UJ-4r-i2w    •.•.•.•.«»•»«•>•»«•.   ••»— K-2:<ril— STHvO— <or>  11     t^HOX)^2:^-'-i<_J^  I   _J 

2      2:           w— .►-<►-«  ►-■►-< —t^—trMO'i  *tmvO>-i2:  •-•.£      >— .0^2:1]  II       Q.O      «-»<  II  OLL       II  •—-—  II  — 

d     wOO<x- — ^-^^^w— ^.www—  aacoaaoozooi:  u  oll  11  numoo«qu. 

Zll     QOQ     »-i<<<<<<<<<<<<<03:u-Q3:u_02QQ     oQi-iO^  Q             Q-^ 
<      2" 

oCO              •— •                                                            O     0s          — *0  O                  <— <c\Jrf|^0 

0.1—2:                                                                   o     >$-        intn  — t                _«     ^     _»        _i     <m 
oocx     2:                                                              -< 

ZKO       O 
LUOOU-       O 


119 


O     mo 

O 

O 

-» 

X 

X 

o 

00 

i/) 

• 

al 

C* 

«»» 

in 

UJ 

UJ 

-5 

•— ( 

> 

> 

•• 

UJ 

a 

O 

_J 

•» 

•— 

B 

UJ 

UJ 

< 

II 

z 

Z 

* 

• 

< 

<r 

o 

•» 

-j 

_i 

1— t 

UJ<— < 

a 

Q. 

»-» 

h- 

Z5-. 

<\i_l 

< 

_J  »* 

_j 

-J 

rO  »• 

DC 

<-« 

< 

< 

••-J 

1 

>— ' 

o 

O 

>*  — 

*-» 

•»  •• 

1—4 

t— • 

ZcO< 

"5 

— . 

y  m 

t- 

h- 

2  .»v* 

#. 

*: 

s:u_ 

c£ 

a: 

►rsj— : 

z^ 

z*: 

«-o 

UJ 

UJUJ 

<-»rr>_j:z^ 

z*: 

-HO 

> 

>S 

a—  ». 

-^ 

»•  •* 

■> 

»-< 

-j— ^. 

— *w 

^^ 

—^  •» 

-8-  UJ 

UJh- 

11  _» :*'  cl  < 

n  ^ 

C\J- 

z 

z 

1^    »>^_J  II  UJ^isi 

m 

Q-4 

HH00 

-^^^^ 

II -O*:  — 

•MM 

Oc£ 

Qi> 

+  ^y<t-nz^< 

QvOw 

<UJ 

< 

_j    ^:  it 

»•►-< 

# 

II  -I-Z 

zss 

s:oo 

II  >t  —  Q«,*C>- 

~o 

'-ulu<o: 

►— iCQ>— i 

CDQ 

«-"CO<'— it 

<1^Z>J-Q( 

DIDi-SZJ 

GOIDt— 

3 

Q.      —  h— 

i^O 

i!  : 

3-J>-c£H-a>-</) 

GOO 

joll<o-ood 

ii  <ocouj; 

Z^O      oQ 

2: 

0»-iQiQ< 

QQQ>2U.a:iU    •      > 

< 

H-> 

> 

— i            «M 

co-to-« 

<\l 

OUJX 

uix 

ro         ro 

r^cO-4-vj- 

in 

_JC£^ 

cto 

Qu^DQ. 

Z>>- 

^OiUJ 

c£  ►— • 

V.QQ 

oa. 

in 


120 


oo  r- 

<m  o 

IU  o  in    o 

Q  Xc/>  O     ro  CM     O 

il    oj        i— — •  •     is.        mo      • 

>-t<'-0LU  I      o        oo      i  in 

31      tOH-  m       ||  lAO  O    t       II  rn 

ZLU       I—        —O  Q\       q^CU    •  «   »         >o 

<o_:lu  a  --<  oo  i  H  it  *_o 

o^:_j  a:z  x  mo  aroii  >co  ►  • 

U.oO<t--t  o— •  «                     o  — io  oo— «<nj  >o  h-  II 

oz    i—  u£  s:                  <  r^  i    .  x  ^aco  ^j-vi  o-?: 

UJl/ilZ  <X— I  <  -*UJ  II  •-  II  — 4  o  »  »•  J> 

</>      2Z5  UJO  a.  O  -O  I  -<r.-»  r\ICL>0  fOin  uj  :<_ 

h-O1-*               eg:  tD  h-  •  Uj^.l<  0(V>    •  oo  O   •« 

2T(—  <H-  XI—  O  O  CMT\   *.  OX  ••ll  »0  »m 

lu    oa:  i— ^  a:  <  ii  c\ia>m  •  -roar  |    ,  coc\j 

5:0     <  ao  o.  'ii  o  o        o^:>  ii  cr-oa.  ii  i  »  o 

aiLUiut—  ujo  </>  <  •    orsi— <o  :soo>-  —mi  ii  o 

oChhifl  Q  LU  Lt.oooro»o«o^:o»»'3:ooa.o 

►-»o<           uj  — «x  cos:  o  —4  •mii'-<o~<ii  o  xo'iim  — «o  -<  • 

OLU1— !-       ai  >*l-»  •«  O  O        I  O— «l_>  I    ||     'U.    «        ►oi-a      rsrvl      Q.  I 

O^itsiO     Oh-  -ho  -fro.  •>    lu  il  i  Olljo  ii  a.  i      in  »>-o     rsj  ••    — i  n 

uuuj     <_:     u_  n|z  ma.         o     ocouJu.aMjjcr3a.  u      so  I    -o      -in     ^t— 

□_</>>-  •  O  I  >-h  —«  •      OO OOufMJ.U_OtMvtO  II  s0    *      -4— i        «*o 

a</>  za  »   t  vta:  o    Hu.oo"ioiio<ooQCon  ~«oa>-oo 

uj^yxijj  o  i  _jq  v.0  iioooin  *^uo  ••  *•  I    »c*:  »u_  oooo:> 

^►-•tJjQ  >-«lo  •  <ulu  I—  j.    troo*J**0<\J  -OO— t'J-Jli  X  I  O  i   iO  i  ^ 

•— «      I— O  I— —<  _J2!^  X^O  <IOH    -a\r-t;\iCO    -   »CN>    ••  II  >  I     |OUJ    » 

h-h-<-0  ■<  iuz:—«  o<   •  a.  i|    «ocr>  t    mqoo«oxm>  •»  u  »    •Ou'> 

Z     zr  Ocx  ^<2T  *-*      — '  Q<  I     •    »uj  j    I     •  ii     •COXOf-trf)  -^^  |     iq 

IMLUUJUJQ  C?Q  ZIh  O      a.''.')  I  OcO  II  O  I  I**  II  LUO'il  I    O   -0>0  -vl,^  II  f~  I 

«ISIh  LUOi  <0_J  II        3H  II  -U.O  II    liyN-ii.miMJ.il     •aW-J0  -..-via.  |  uJ 

»—  aiv- h-  a:  X  lu  11      k-  vfOaD<ujr«-Q^<,-L.-^mr«-  «>-o  — *  -ex.  11  o 

2:     a:    <  zijj  oo>-jj  lo  01-0  — iou.uo^u.olio<  ii  o  i    »  •  o^i^ik  • 

<r_j'-'  »_i  — «  lo        <     i-^  o  >~<tLa_  1    ►Qu.u.oo^iO  »»^sjo  h  — <  1  00  -Ovf 

a-»-ooo  <x  3         •>>-«<  <  s:    — 4  uiO'jQo^jma  ~:n— «  •  ar^n  00001 

OH-  3TUJ5;      Oo      ^>  I— h-O  •>      ►—  ?.Q.    (M   *   -'^)L}    •   ►   >Uf<l^   "I  -30I-         •    •    *CC\\ 

QZujih         h-     _j         u.«Oa:^)     — «     -JO  1      000  •»  *  1  o  1    *».-« <*•»!■  li  x  »o      l    lO^a: 

QiOcjiO'^      LU»-h       _J  ►UJOS       il         I  oc  1        m  II     •'NJO  II     •  II  O  I  NOOI  •»  I  O       II    II   II     «-i£ 

o.        <        ia  o  j-oluoto  h-  1  a.  1     n.j.oh  •1qoo  •luoj  i  ^.oho.  00^2—4  ► 

Oluluijj     h-  u.  <^>a:_i  <  1      Q-     — <<  11   1  o-O  ti  ujo— <oujx(nj>>-  N>:<or\ 

LUSICIX           cC  I  iJ-»      LU-J  h-  C~o:£      tMij_OiJ  II  ULUji-U  11  TOh-   »»0>   *  »:siNOO 

»—  >■-*     h-      ^O  i>1  I  >— I     O  u_  ujrxo       •OaaroQOOQU.m   •  0-1O  •»0">  f*  *  ""O  « 

<H-'-')  •-«u_  <  I        <XU-  •  O.'Ma.O  il  OU.-4-OUiJ.OilJ  O  il  :3-^    «^\0  QvOO    »UJ 

_j    — <  «o  oo-j2!t—         io"s:o      •<  »;3-n'i.  -o  ••u_-^_i  00  i^o     00  •  1  s\ 

O'-u    o:uj:^k     o        j<ioaw      1  11  oa.u_— kcunocooo^'S  •  •1100      iqoii    • 

OIZU>OI      LU  LUO>~tUJ<      I— h-O.      Q  11  LL    •  ••— O    •    t.  «o    •    «-in  |<  .-ooo        I     •  II  TtlfO 

_il— O— lUJi— O      Zl-      >-«l— O  <S.<t      O      'UOO^    •   i-O— «0   »l  roll  C£0    •    •       II    I  — •  J  ♦ 

<  — »_j>--<o>~<    oa£    iui— <    o        h-ai  «2!:  »o  11  — 1 11  %*■  11  —*  11  m  11  m>->co  11  ii     Qlii"^^;ii 

O    »K-OX  ^UJ       —•— I      Ja<t-^Lli^Q'i.fO^CJO    -U,  1   a.— <<I  I  CO— iu    iwXQCLH.       >— ><-*    »-a. 

i-'-'jJou.j     loo.         ujo"^i—  11  zz  »     ii  — «  •-0-tLu.xi  1  Qlulu  1  inn  ••  -x  ."i  D  .->j  ra^^o  ^: 

ujujoi— <            lo  1      ^>-ici<a.<— <ot— t— o— 4'J-(\j'-Juaja.— «;j.iuLLoj  X)u  ,  »->o  »— <  *-o  •• 

aiX^^      Olu  •<  1       2             -^5:      11  t—  OJ>o  1    juio^OoomO  11  o— <in  *.j£  in^fNi  1  ^ 

<I      OOi/)ZIu       I  Q^<JJLUuCLUOi— '-u      .tiLuO'-JO— <'->-UO— «0  OOo  I  O'NJ— ':-  n-^joluO 

(JUj^-^u^h.     ujiucli—  IXoaau,i/)k  — 'co-o-n  ««rf>  #»o  •-^r  ••m  »dujoo^  »>  o  -ooo 

CQI-U  I-       XoTi—  >ii      H-H-a'O  I        0'-^j«-«  'I  OOinO— iO(^OChO    ttOOONiAO^O    «0 

s    z^ujimi  •>-<<  1  s        c£      1  zo  1  odj>  m  «-*  im  »h  »oh  •  •  «oc  •+  ,;^  • 

)     ,»«<t>-f_l— cajJ-*      .£  1  ~«<2;zaji/)  I  uj  |   3  I         »XO   »0   »0   «0   «0    t  -»  |    |    |  OO  |  o  I    I    t 

<  <ro.<*>~H     ^S'^>      1  ii— ■— *     — «ci.lu  1  m  1  ojo^  11  11  11  11  11  11  H  11  11  -.nm  ;i  11  11  o  •  11  o  11  11  11 

•       Ol—  h-UJI—  lUX        0C.UJUJI—  X        S"2    I    <^l^<     •     *OCuLUOU_0<UJO."4     •H--/)l—     •   I    V~     •!—►-> 

>  Z^O^iJQ'UZI-       U-.UOO.Ol-X       1^0lK0  3300Ci.:JdJU.^NOQ^I    HDIDQ-t 
.       OLU  J-LU^-  <<  —1^:2.1—  "S^-JJ IO"S:— i.U— •-ju_ii.LLU.LLO.ll.U_LL  II    K  Q/)Q  II  >Gll  QO> 

:  <;2:>-LU>.Luajt—  .jj_jot:-aoo— 'O     a3<oi— >  11  .ujD3JjGOqoo/Oo>ij:a334-i 

<  saiuicisxza     -•  >i<^ic_a.f-        siotuoooooooooooo— «xx>>>^J^J'^>ii«c 

»  ,-4  h.  <  i—  <  i—  <  o  -3     a.H-00         ^aaMiiia. 

j  *~<dj    lu    ct:— "•jix-Z'-i  a<ujiu<t.5:5:x.i:     zs^zzxzsismszkzs^xssss?; 

>  Qui—  I—  Zl-      <iu<OZ!iiiUTX^DD<''-<CM<<<<<<«<<<<I<«<t<<<<4< 
j     Ooouo^^     o.-^iZa.'-'QOi-i-^-a.ZiQHCL'CLCCaiXiXixy.jjiaCiXa: ^^ctiuctxa:  x 

.#  #  **  *  *  ^<nu  iHf  -^  ^ft  #  »#  •>  *  *  #  *  aa^aQ-i.aoiQ.Q.iiiaQ.aaCLaaaa.3. 


n 


121 


in 


o 
I 


II 

> 
> 

•o 

mo 
oo 
o  • 

•  I 

I  II 

III- 

^o 

~«0 
3:2 

-<s 
s:  * 

••IP. 
foir. 
00 
oct> 

•orMiJ 

II  «oo 

3:  1  o    • 

2:  II  Oc 
~CC    •  I 

mnc  11  » 

oiai- 

00  »-a. 

•  002 
>$■»••» 
11  1  o^ 
</)  11  no 

•O.SUJ 

nJ-s.a3:~«o 
O'-t  ^s   • 

luoo  n  1  -J- 
mo^o  11  -m 

•  •  «000 
-\1  II  IOQ.  • 
II  Q.  II     »2  M 

a  or  o  h   ►> 

••-}■'-«  "sr  i  «■■• 

00  ••mo  •> 

I  olo.-m  .ir\ 
uj  «c\;onr- 
O  I  OO  I  O 

•  II  O  *  II  o 
h-f-    <  1  t—    • 

II  O  I    II  O  I 
dan  i/o o  11 

<  <  <  <<  <i 

cc  cc  al  cc  lX  c£ 
<<<<<  <r 


O 
O 

• 

I 
II 
X 
a 

•• 
m 

o 
o 
o 

* 

II 

I— 

a 


o 

o 

I 

uu 
o 

• 
CM 
J 

II 

Q. 

<i. 
»> 

O^Ni 

00 

*z> 

t  o 

!l      • 

ce.  I 

-+:i 

H3 


o-< 

00 
•   • 

I  11 

II  > 


m 
o 
4 

UJ 

— <       vO 

v^fNIO 
-«C\IJ.J 

O— «:0 

•   •    « 

-*  11  ~< 

II  H.  II 
MJT. 


o 

1 

LU 
OO 

•  • 

II  II 

H-O 

a< 

Q  ► 

a.00 
•«  1    • 

OLUO 
•O  II 
O  »rX 
II  (MO 
I —  tl  — - 
OHO 
0<    • 

O  »o 

oc  m  11 

O  I  Q 

•  UJ    •• 
OOOvt 

I!  •  «0 
KNO  < 
13  I  II  II 
O  II  CO*-* 
<IH-  O'vJ 
>U.       ^i 


<<r  -t  <t  <i  o  -a  n  < 

a.  a.  Q- a.  3  .-<>-«  •-,  o. 


UJ 

00 
O 

1-^ 

co<o 
ot— a> 

O'-L  •> 
o  «-oo 

•  CL_) 
II  TTQ. 

a 

SUJUJ 
<5  O  i3 
(XN-H- 


ST 

s: 

»— 

CO 

O 

a. 
5: 

>— 
00 
a 

■£:< 
SO 

aoo 

00. 

Q.< 
S'— 

I— U_ 
COO 

Ql 

X«X 
lUQJ 

00 

I-  I— 


if 

CO 

_! 
LU 

Q 


OOOO 
•J-J-l-J 

>•' s^ 

^<J  <Mc\ir\J(\l 

OO  OOOO 

^^  OOOO 

l-D3D3»>> 

o>*  •#  #  ■»•  zzzz 

0>»>WHMM 

3:^:2:2:^^  ><•  ■>(■  * 

21  **h  •—«•—(  >— 1  *••  »^  >-^  --^ 
*  *  *  K  ihU<^3<^>^ 


O  OOO 

•  QJJ 
II  C»^^ 

y*  f>i3:3 

-j  1  -cs: 

UJ  UJ  #  -^ 

Q  UJOO 

•»  <lUJo 

0  -r-  <t<r 

•  hoov-u+  + 
o  ojjoj?^ 

LO  o                Q^-tt  Q-^J^ivJ 
i-O  (MfM— I XX  >  2  I- 1-  3^  -«•  -fi- 
ll O*  *  tt  «-<M.M002NUJtU 
-h^:  ^ ■*  -M-  OvJ  I     1     I    \    QQ  I  O'ilU 
K-  II  >-<  0000^>rvl_|  3:C3>-    •<< 
Ql-K  «0— J—1  m*i-«»-i  >"  >„ -vl  >-«  ^-<  ^-^*- 
jOZ  —>  II    II    II    II   II   II    II    II   H    II   II    U   II 
a.-)-*  II  AjrNj,--lXX>-OLUOUJ>^C\) 
^Q-U.  HjOUNi>NOJLUUZ^H 

OOOO  -J  —•>-"••-•  <i<<^  -i— <<< 


00"^v.^"v 

00  '-roocri  co 

WQOQQ 

cycssss: 

•£.T.<<-  *  ■»  ■* 

•«■  -K-  OOOO 

OOUJO  JJO 

-UO<<<< 

<t<  +  +  •*  -f 

•4-  +  LT.u^CSJ^Ci 

C?3tOOCjO 

rvjrvj^r-sJrv)rNl 

■«■  ■«•  -;;-  ■«■■!{■«•  0000 

UJLUJJUJUJUJOOOO 

LUOUJOUJ  J^^-^n^ 

<<<<<  <— ««-«f\ICVj 

v^<o«>~ »^^— »•— »~ir\j<— <i"\j 

il    II    II    '!    il    II   <<<< 

-^<"\I.~<(>J,— »C\J  II    ii     l|    It 

(\j  .ni  ^^  .h  ro  r\j  — <  r>j  ro  ^f 

<<CQCQcQ'X)'XllQaDca 

uji— o; 
o<^oi— 

HZZ 

ZQ^ 


_J< 


?r 


U 


122 


< 

+ 


oooo 

^"^^^ 

-*_*f\l(\J 

(OcQUOcO 
II   II    (I    II 

'3D  CO  CO  03 


DO 

I  I 

II  II  ui 
<Mf\J> 
>-  -4  — 

OOH- 

>a: 

—  o 

OioO 

QZ 


Q-.  OOH-  «^.«.~_ 

Oh-  OmO  —<r\|rO-$- 

rsJO  QtXQ-^ -•'-*'-«  »-«--< 

oibf— ki—  *—>-  »o>-qo^:v:^^ 

•  "QOOOOO    •   -GJ3.   ►  »  -  ► 
Oi'WOOQQ    «00   -  -OO 

— •i*v>:3.o.  3»l^o«--  »:)o  «j 

__j    •    •»   •■   »  v  ►<  jO    «    •  O  — 
C^tA    •    •    •    •    •  _j:_£-~.OQ— *-«»^-'--» 

OCNJOOOOO  HJJ^-JJJJ 

if  —I  -  J  .  J  -J  _J  _J  H-  <2  O  ,-i  LX  ,J>  O  O  o 

—  a;  a:  X  &  ol  vl  ~£  >->*-  'S>  S>*~  t- 
ii  ocsooo  .*)  •■*  ii  zi— i—  -zsr 

I»-KH"  I —  ►— ! —  i|  ."£>-'  ^2T •"-"-• 

KZZZZZZJJII«-'II   11   It    II 

Q_«— II— IH«MM-<  _Jh—  -J   '  I     II    —*  Ml'O-J- 
LU    II      II      I)      l|      l|      II     0»-H<£  — (Hrt-Hr-lrH 

QD>^a  Trac'-xavu  ^^ii^:^ 


r\JrvljT| 

aoo 


-3 


ooo 
ooo 


OO*: 

»•  »-cvi 
OO— " 

•  «c0 
OO-*- 

— J i  - — ! 

OcJHJ- 


^"  *  *:  jl 

**  -ir  •»■   £• 
.•\JC\3C\JAI 

OfOt^Jj 

4  4  4  4 

c\ic\ic\irn 
<-  i*  -ft  » 


^^^ZZ 


II   II    II 
txjr\lro 

^^*: 


^^— «f\J— trM-* 
— <  ►-*  O  CO  CO  CO  o 
II  II  II  II  41  II  II 
nT  4"-«'\l-*<\l-* 
r<")^*-*--«r\ic\im 
:*£i<;xu_u.u_u. 


«■  -»r  -If 
r\)i\lr\j 

cGdjro 
4  4  + 

•tt  -«•  * 

-*~<.~ « 

JO  ID  CO 
II    il    II 

U.U.IJU 


c\J 
«-« 
u. 

CM 
CM 

Li. 

+ 


-»  I    CM 

r\Jc\i— < 

cm*  rvj 

X-4'.L 

*  \£  ■*• 
<-xj  .*.  . — 1 

■;<■>-«•—« 

■ii-  — <u_ 

— !<■& 
»-<■}$■  1— I 

x  '  vj  (  n 

i't  ^^ 

-*+■■& 
O— 1 w— 1 

I  ~"0 
HJji   I 

II  II    l| 
r-tc\|.r> 

OQO 


(M 

# 

<M 

c\j 
U. 
+ 
(M 
■K- 

(NI—<(\J 
•-•OJCNJ 

U.LLU- 
*  —  # 

vt— >m 

LL'JU. 
+    I    + 

r-l.-\)r-t 

sjTNlrri 

*o* 
1— 1  ••— < 

I  +    I 

#-*  ■* 

-tMv}' 

+■#  + 

COO-M 
•-<  ».«J 
^M<[ 
+  4.  h 
r\Jr  "vl  r<-, 

<:*£*: 

-!«•  #  ■*■ 
.— <   ».-« 

II  11  11 

aoo 


>— < 
o 
J 

CM 

CM 

+ 


CnJ 

< 

CnJ 
N 


CM 

u. 
<\j 

I4_ 


r0 
4 


CNJ       -*( 


OO 
I 
AJ 
CN/OJ 
■K-  <T. 
■&  -> 

— <-v|fNI  '-n 
"• <  ^J  ->{•  CI 
<!_>«■  i£ 
■«•  «•  — I  + 

:\jvra.r\j 

4-  +#  ■«■ 
— •(NJO'N 

*:u-  »  ^ 
n  -;(■  H  11 

<->j  sj-  r*1  .r> 
OX  3Q 


1      ' 


123 


CO 

■f 


•4- 

LL 
4. 
(\J 

* 


» 

t\i 
t\l 

< 

•4- 

O 

4. 
•4-0 

x: 

*o 

ot~ 
•o 

<\JO 


MO 
<    • 

O    t 

•  Ui 

r\ix: 

ti  >-• 

Qli. 


o  oo 

•oooooo  •  • 

&»•••!      lOO 

H  OOOOOO  II    II 
»—  II   II    I!   I!    II    II  V-l— 

OHhl-hl-l-JJ 

QOOOOOQQQ 
OOQQQQQO1- < 


e 

00 

UJ 

z 

< 
-1 

Q. 

3!3 

oc 

aa 

UJ 

QQ 

3 

V- 

*»♦>—» 

OO 

OlO. 

o 

~<f\J 

0 

z 

•*>t 

o 

< 

>-> 

+  4. 

oca 

UJ 

CXQC 

< 

UJUU 

aa. 
— ir\j 

»~« 

>* 

rOrO 

< 

>«. 

>-> 

u. 

---to 

+  4- 

LU 

|  --4 (MOV, >-•"»-» 

-4— 4  GO 

ujuouuuui- 

II     •    1 

nC  il   il 

#tf  tf  #.*  # 

* 

•,*■ 

•ZcCcC 

II  OOcocD 

<—\ 

^.■-»-«» -»—«—». 

-■ 

--»iX 

uU.TitUr^OaOO 

OO 

fvUNJtrytrxOxC 

1^-S-O 

O-'ajOLUr^l^J    »•   »• 

^->-o 

2SSX2S 

^•1 

STUL 

CL^JO*    *   r~i^ 

OOv- 

rO      <*>       rO.->(^Jr>J^Nl~J.vi 

-,' 

>>l 

-* 

M-^fNj    •    •££ 

ooo 

s\     la     en  •:>■■«■  a  -;<•  -:'<■  -sf 

-k- 

*aC 

.t—-. 

rsio^j  ••  ►O 

o 

o>  0>  o  <J>  O  o>  J*  r»  co  r*»  -cot--  co  r- 

CO i  JJ  O  i£  1*-  <*•  O 

-^>->OOI- 

-•*t*^« 

OO    «,0    «0    luD-JiQcOlDaj 

rc. 

^j-jiXQ-^h-    • 

o  ->»  *    •    •< 

oo-« 

\O'9'-«~0'-4>O'-«i*-+'4"-^4-  + 

■^ 

+  -JO 

'3      t-IOJ 

•  .ttooi 

•    «-H 

•   »m  «tr\   *ir\~4~ «ro  -o<r  t 

into  io.Nl 

rNJ'-iJ  j  («-»-—  OJOOcQOO    • 

^^r~«!M^-«-\jr-iiNj.'^;.X'5;^:i;;£; 

a: 

>ZX  1 

1         •      •       • 

■OOJJZQ 

Q    •    «3 

-4*^  O^^^OM^-JfvJ-V^ 

N'-^^X 

Xl—  hl-hN  VIO.Q.UJ  i 

1    CJ»--UJ 

o        -o-O-^'Ch  o-*<-  h  *  x  -:»  j- 

#  ZUHJ'JJ 

c#  •»<•  -j  jooaau  • 

t      UJ  ■«•£-#  *    ft  ^  *  ui  O'^  •Or.OiA  OQI-O.    •    .    . 

«^(M<<      QQ    •    »0_ 

O— 0(N-^<\l<n<3-  ■t^OQ.'.U  £)-?3CO-OcOcO<J,,-« 

Luac^rac: 

c>;— i— iluuiluoo 

ajc^crrso— i 

II  <<")  ^ •-«<"-•  I'M rj  Nl1^*!  •*}•—•>»»»_.  ^^^#--* 

w 

—      iOa^aJiJ»--(;^JG 

OLUiJJ—   II     II 

H»       —'ii^^^^iii^l     I     |     I     1     | 

1 

1  O  II 

il  uju)  ID 

L2  II   n   11   '1  <~2  Ii 

II  jDc&DSIO. 

OOI-  ii   ii  ii  ii  ii   il  n   li  ii  II   ii  H  H 

II 

II  »-t.X 

■■x  a  a.  --j  (»si  c;  o  o  o  2  ^  <x  a  o  o  5."  x. 

:.):—  ^T^1^  vvi"in  01s-  — «  »j— «>\i-M\j<~« 

.%;:—  cH 

LU-.^h»  •*■ 

«•<!  ^Ow^'JUjw.*-..!-!- 

<i  u  d  s  s:  s:  s:  i:  s:  s. -^ -<(nj -ni  oi  ,m 

-J- >T  <  U-1 0  U.  U.  U.  LL  "VCD.  V)  03       (/} 

CQUULLLLCO  13 

>-OO-^l^0f>J>«ir>l>lr-J>>>->>>->->'>l«i.'X'NJf^'-l'— < 

-aaoQiuo 

QHWhOJ 

O 

z 

O 

H- 

/v-_^  < 

Itl 

O 
< 

a 

124 


r 

D 


o 


CT»  CvJiMCO 

09 

■  o  h-        a  H- 

OQ!-  OC  H-  ^ 

■ft 

OO-**  -»  -*•  **» 

oo  •  i/>         •  s: 

•  «Qf  O         or  m 

Oth-uj  I  <t     LU  I— <? 

Lti_J    •  UJOO  •  I    00 

cd»-i»—  o-«i—  ii  ii  a^ou" •  ii  ii  Ztuujt—  ii  ii  vo— «t— o<vji—  ii  u  </>>-« 3o<m  n  ii  z *- s: t> cc a. ♦—  i—  rs  t:  u  'i 

ii  ii  iiQoaSiii  n  iiaazzh-QOJi^ti  h  n  JOQii:  11  11  iiQor^hDa.ooan  ii  111— j  — 

U^-<  II  K-  —  hv—  U>-  <  l »  » —  f — » —  2>*— *^i-UI-<  il  h— »HHUh-<  II  i—  I— ^-  it  kjr  —  * —  —  <-H(-><  ii  a.  s: 

cQcoajn  ^llcoodx)  act)  i)(_K;OLO:3u_i_m_  vko-o^ooo  ^LJu.oooo^ry)^^Oi/r-oouJ:_i.'~K;_u_  .su__>  j:3^> 
UQQQOHOUaaOOOQOOHMHQOaQOOjHQQOQQQOQQu      —*'^i— «~.ii.  <  <-t  ^  O.O. 

Ci. 

oq  h*  O*  CM     —*  ao  f-  <^  r\i     — <S  «H 

t-*t-H  — (  r-l  -t  f\Jf\l  (\J(\)  f\i^> 

a 


-J 

O 

a. 

o 

# 

•^ 

*^« 

i-  . 

.-* 

CO 

• 

a 

Or 

1  <s 

UJ 

LUCO 

* 

LUSICj 

s 

.* 

_l 

o. 

# 

~> 

LU 

s: 

>— « 

►-< 

1  D 

uji—  a 

a* 

fsl 

h» 

a 

cc 

t- 

_J 

a 

TX 

cj> 

* 

rO 

o 

^rvi 

•— t 

O 

OO 

a 

>- 

h- 

-»a 

a 

OCD 

'j> 

• 

t-«-» 

— « 

«  -" 

— « 

•    » 

• 

Or  LI 

i/>c? 

LUUJ 

ZUJ 

•    • 

UJD    • 

H-t/1 

O  — 'X 

O-J 

125 


o 

o 

h- 

4- 

OCO 

o 

o 

c£      O 

f- 

0."»f- 

o 

l^OO 

o 

#~«0 

o 

—51 

•«»  • 

OS  ■-•-*• 

*-«  — 4 

•3:-J— < 

•  + 

— «  |     •    • 

or  2: 

■*ixk~uj 

UJQ 

o-sioz: 

•  z> 

:di-i  •  • 

z:o 

f  jj_  -~<-— 

— «o 

o 

a 

O00  o 


-J          X 

-J         u 

O          h- 

aC         kh 

w        a 

^        — 

•-«        </> 

i/>       a 

,*             !J 

— .            \ 

X          -» 

o        — 

K-      -»_J 

-1       -J-J 

a.     -ia 

«-     Oa: 

00      QCw 

o^ 

o    — 2>» 

LO 

o— z^x 

.-* 

*    JMW(J 

^j- 

>_J>0-H-t- 

~4 

+  a  -*  a"-« 

• 

—  atat +  a 

m 

X~»  |  -»»•» 

>.. 

o^-^-jz 

X 

y-o^jw— 

o 

O'JU 


'Ji'-OCO 

or     — Q  O  X"-*  I •*  v«QQDO  rgfMoj 

C     Oh-  t-  o^-^-jz         I^*-»*OOOOOJ 

s£     ►-•a  O  h-0_j_j—<         or\ia:ooca_J_J_J-j-j_) 

—     »-'  Za-M-  CZcCii  oC7>'—r---K-  *•  •«■  □£  v>  dj  or  s/>  Q 

ix!     _i—  <-•  o-~"- -cc--*!— fvrvjairoxoaacjinocs 

T5T       •>-*  «-*  h2I«*(/)Q  «  i*#  -::-  *  #  ->  -:(•  ■#  -ss-  #  # 

-**-*UJ    ••—  — •        •  •-«  J-hO^OONN    •-jO."D^Z)t53Z)  J>0 

h-xoujxx    !ux  < 00 hQo<tn">ojit  ■»■■**  «■  #  ■*  ■»  ■»$• 

ijuo^a  1  •— 'i—  \  +     uj  +  oo"T>i-U  1  -'^aOQii  11  -x.  '-n  ~o  ■£  <.o  cc  at  1/1  ao 

o^»>:<rr*-2:  11  >—  f-"Ds:  11  h  —  u  x<ri-^cM^X7>ooo  11  11  z  11  w  n  h  11  it  11  «odoacaoon 
11  11  <i—  r>.t:u.<<3:"-'!^^:  n  i.-i>- j:u-  33  u  w  » ^n^. ^vt—  01—  n-r— << at  txx;x>  si --jest's: 
Q.M£itOQ.r3il  11  iiH.jw«»-.,i:u.u.<(-"D<<aai-Lii.j-"-a;Qij:'3GG3:a:ojii  n  11  11  11  11  11  11  11 

_)X--l-«-^i-KX  II  iU.51X><  —  •—  II    ll  I-— »  )i  1—  S < '>") (X  >_* Z <t CJ ■*<•  OQQO  0»--.J»-«i"\l  '<">  — «c\j  .-*■>.»«  (Men 

►-« 
r\i  ro  >r  Is-         00         >0  <7*X 

Z> 

< 


126 


og 
o  o 

-J  -j 

>—  -*  • 

•  #  CO  CO    • 

«  Q.r*j— »  S*"*  • 

•  *  -HO  O-H-X  • 

—  X  2£Itf  *  ^#H*  * 
15               <_)                I  DK-                        — .fl>>-<                   -*» 

<  -»(->  3  +  lL  -*r>Q.  H- 
+  ac*~>  -K-3+  X  +  -—  < 
h-  «•  a.  x>u3»~  o3;v)  * 
u.  r>^>  *  <r<t  h->o  m 
+  »■  to  -j>->-*  *-^<o  * 

I—  -JO  S— •*  Q.-B-*  <y» 

<  2:0  +  3:-»     u  — «I3-<.  .~» 

w  I  *  O.r-1-J      J  2:-*H-  t\| 

—  -x- _j  -o#o  -*^J  to:*-*  h. 
X        3:_j           <Z3.i'     -~»o        -it  ~<sm              u. 

"»0  -K-  O  "r  +—      X-J  CD'S    •  # 

-»C?K  _Jct  2S'/1      0*n.  rvi  4>  ^  u"\ 

*  3:0.0  +z  <<o-*  — <^m  »  *-<m  m 
oi*-»^  a>—  #  •«■■«■  c?a_tt  a"!i  r- 

*  o.?:-»  >to  c?3:~»-*  •—  ^  r>  #  >H»  — 1 
a.  "v  — •  (M  cc  *  -th:-  x  >-  to  :d  -><-  2:  <  u.  — 
ac  1  to#  <o  <-«jMa#T«  >-!!•#  o     + 

+  8  00    3:^  _l-^'-<|  #1^0     -fr  — 1  •        ^fM-rt- 

-«.-.  '03:03      >^  02Q.->-^  +  J       TC^m  —OD 

ooiH:-\ii     o    3    — <m— a.  o  3:  -^ac  — « —*,**—        ar_j^<M 

•  •#  :?--». n    * -*     *-     s  +  y>#  -j-K-  >  *  ^s:-h        #>v>o 
~h— «X3:rvj<     a;^s     -»     r>d<\ioa-  Q>-tf  ~»-<  +  --»         3>2:o 

-  ••3J- x  «■  +  32 >  o— +  •x-o.-Qct:^  >q_s:m  1  >•*-—  v. 

aoa:  +  *o  <<  _iaa:  *  +  z-~  -s-^*  -*■  *  -^  <t^.v-» 

•  «a:if^3  & -x-  ^-x-  *•  >o-v; rZ>  +•  x>#  *  *-!o-. 
00x3- -»■<  — «>  >:x>-is-  -j<-v_v:o-s.'  :•  3*3  a'>-j— 1 

-J  •  •<}  cix^i    >-<    >  +  >si-*-«.c>-.or^  ^s    «■  — «>»■  oitototo'^; 

^»  to                              -~i~*X  >  +  <C?>>      -t  CiC ->i  +  3^.^j<  j/^<.  CL^Z^cNia.  3--=-%~-+ +  +  + 

<Xro  2:                     —1  1  cyxtN  ■*.  «•  +  4.     ql-(<.  +  ,-nj><.  +  siic  >  ^>  •»■  +  ow-  +>orz:  ^-h<^ -n 

o<  o                     3:  *•••«•-♦•*  r<.io.ct>     is-  ^a.-?f  ifx  1  >^_YtMoJ-  o».~<<-  >—  4000  c 

•«•  a.  »-h                   3r^3>  ^n^^  ^  ■«-  -f  >     a.  ji<- «..vj-{;  ^.{'c  ^ ?:■*■* -j^x-i:- zi.-^^tT^^T^ 

Z>+  K                    ■»■  >—- •or-*  >-s  c?2.4;— >^>  <+  J-;'-  >+  +  ^.^  ^-- >2  +  "•*♦  +  +  +. 

■»  r>(  i-o     n     <t-*-s'«.-»-«3>3:^:x">-«-  n>a.-«-  -n-  Z)*.-—  *  -js-  (j* .  s^aia  •:<•  >c-»—t*  "s  +  — ..N'rxjrvi-Nj 

x»<  o     uj     _j->sglC?3;  -f-^oooo^--«-  >^  >ax^^  a.^-;;-  ->i-;t  >•;:-•»>  ct:*  «•  ->  l  "±-i  aX<muo 

*a.  o.    a.    uj^>-^--»-->h-^7: ■){■_)>"<■»■  ^?  2 +  v  cr:>^o -?a.'-^QLa,>>  D.J3  3> ^  oco-t^ijo 

of  +  4.     +     oito^ooooLo^-  xoJosx<to»"^i;«  ">ifNHf  v:  5^^:  ^s:>^  oaz:jau  +  +  *  •*. 

2:<ua<_)Ovuu.nr.  <<<i<  <ny)ii  ii4.  +  -*-4.  +  +  +-f^  j  +  +m  +  ^.  +  ^->-w  —  2i^»-*  +-»Q.<rjOQ 
it  o.ao-aaa<  11  11  11  11  11  ~s.  11  cxq.  11  h  u  u  11  u  11  J—  ii  11  11  11  11  11  11  11  11  11  11  11  11  '/>  11  ;,  si  11  w/).o-ooo 
— »  W  ti  II  II  11  ll  iu^^a.^a:^^^>^c\Jr^ -f  ^r\j(^<.r^^;,0-4-'-<Air"i  4-^-«f\j -^  4-'_j^Hf\j^-)-t^o  11  il  H  II 
U.<^UQUU.^I]'QQjJa}>>a  3<<<l<l;0-0  "QwO  OOOOO  QO  13  IUUJLUUJOU.O.U.  UUQ<t  t30C 

txa.o.a.a.a.u.-rt<<<<o<a<<i(/)i/itO(.oco<.o^o^  oii/moto'-O'/j^'^otooTLOtoi/x,-  (•)vjoTwO«_;<.o'y)cnto 
-J 

2: 

a 


127 


LULUUJLUUJUJ 

oooaao 

U.U-LLU.O.U. 

****** 

O.U.U.ULU.U. 

U.U-U.LLU-LL 

oooooo 

+  +  +  4.  •*•  + 

HJUJOJOJUJOJ 

*  *  *  *  *  * 
UJLUUJ  iJJUJLU 
<CQOQUJU» 
U-U.LU1.U-LL 

aooooa 

OOJJ'JiJ      > 

+++f++     > 

QQQ3QQ      '-U 

******        + 
OOCOQQ      <NJN 
<CDOQli-'i_V       *  * 

u_u_u.u.llu_  ** 

OOOOOO  OO 

OOUJUU  ** 

+  +  +  +  -V  +  -»— . 

OCU  JOJ      cQCQ 

******       *  * 
OOOOOO       O'-O 

LLLLUU-U-U-        ■♦■  + 

OOOOOO     ^i^         -»-»-.  «. 

++++++     rc-M-         ujccot-o  ujo 

ca.-occ  races;     mt^        QCiuuj-^-  c£u~v 

fx.'."-^,!-^^^      lC.-C  a.>>      k\i  (X      r\J 

*  *  *  *  *  *     — '  —>*        -  »•  -a-t  oQin 

CtJcOCDrOrOca  +  +  OO^I-  -MV— 

<ffi'JQ'JJii,  OQ  •    »MOO  >00 

sf^                                XL.U-U.'i.UL  «■*  OOQOr  ac                   oot- 

UJU.                                 Q3C20a  —  ---  rqrr>  ♦■      O  a.                     •■     O 

</>«/>                                UUOOOO  — «•-•  ■*  *»f\f-»(_9  — »                  D                  c\J*^0              -»«^ 

*+  z+  +  +  *  +  +     30       ooni"d  a:  a.  >-o  oo'^o 

raro  0<<t<t<r<<  x>*  ■>  •   »o   •-«.  a.  <  Q   •«,  O.Q. 

UJIL  *-<^>>{.-vi'-j-J.-J..jr-j  J-  OO^O— »  3C  v»  --» o^  2"5I 

<^>t/>  h-****  **  Otcco        w—»a.  •«  z>  oj  a.  •  •  03 

+  +  o<i<r<t<<<^>.+ +        jjooja  o.  a.  oooc?  a.a. 

'M'NJ  S<^OQUU.>WH  G-aQUJUJ  »-  s  QLUUJ  -*■— 

uju_  u.  u.  a.  ijl  u.  u.  >  s  3: —<  -J-J-J<T    ••  (-00  ID  <♦•  OV-O 

W^<.ooQiJJU.U.Oar30CJQ<->ft53,«'-U^X  cq        •  a.  UJUJ--><£  co.*n 

+  +Q,(vaaaiQoouolJoLJ- n>%  1  )'Uujqido-ioo<     —•         1  -<o     DQahnoo« 
»■«—«♦  4-  +•  +  +  +     s-<--v^i3a-i3a:ti:  11  s<  11  11  11  11  u  -« 11  -*o  11  in  n  o^  11  sru.  11  11  11  n  11  — « 

UJUI_<  OU^UJ'XTI  !|   11   II   11    II   11    U  — -- — -  JI   II    II   II  QC"D  0<X  -IZt  II  Q  II   » »  I —      X  II  •-«;/>  30  ^i/)Z'/V-/>  II 

oo^t/ioooooooot/^ziKi—  ►—  i—  s—  h-  >-n  1 1 1 —  i —  ^0:0.^0.1—  z>  jHa.a.MU.K-'jojai-  aaai-jnaaa. 

ii  ii  ii  ii  ii  ii  ii  ii  joaaQOo>"'-ta:^ijja'£-'--u.<o.>:3ss:^<n:K<iszz-^u.asD32 

OIU.<IDOaxu.-<QQQ3QQiUQOaj^OC!iJ2ij,u.DDD3aDMjHH333a3U.a.JDDaaD 

< 

Z3  00  O       —»  CM 

LU 

* 


128 


o       o       o       o 

•  •  •  • 

h-        r-        r-        h- 


m        in        m        in        m 


QC— «. — .(X^o  «J  _J          _juj  -J  _| 

c<ru,<i:o^a.  h-  t—        i — « —  i — —  y- 

•.•.«•.»....  2:uj  -z.-z-z.-z. 

CO                                       LUUJUJLUIUUJUJ  os  o           Ooo  Q  o 

\                                .-ir-tcxjrn^-*^-.— (     cl  x">  ix         0£O  o£  a: 

cO                                     •>•«•>».«.*•*.     z  y—  *—          i—  z  \—  h» 

!E               mrOro          ^v,..rf>,-.,«-..      "x  LU-J  UJ      O'JJ  LU  liJ 

O          15      O      U<*)U>          OOO-OOOO      <X  OZ'J  2        «2UJ  2  Z 

•        a      •lu_jj_juj     3: 3: 3:  :s  3: 3:  :s     o  q-h^:  >-•     o~<z  —•  «-«uj 

~*         \  hooj;-  •»■■»■  3     <x  x- -x x x x  .x     2r  x<  ar     Oa;<  ■x'uj  x"s: 

II  r^O  II  -T>  II  2^r)t—  Z      OQOOQQQ      aJ  ST<  <      h-\<T-_|  <  S  <J^h 

Q  II   III-'      QW<<11^  i-i5"UJ  SIi/JI^-SZO.  3I»-h  -SII— 

»-«i—  luoo— •!—  11  n  hh-     _j.j-j-j_j_j_j     _j  00:0:2:  qd_j     n  coi—  t3 

£<v^acr-5;-2:i-.i-i-^     -j„j-j_j-j--j_j     _j  yo^  r>tu     ^>x  ^  =>^ 

•~<l—  ZJ'-'O— «t_JI~0»-0      <.<I<l<<X<<— 1<  </)CO,J  00>O-^uj  or/1  o0> 

JUKOOJU'Ka.o     UOOOUOO<U  •     a.         aj   •     h-         > 

LU                                  Z  i~  >•  >_JO><  >  >- 

O-*                 '-t-J                                      •-•  ULU2-  UJ       OUJi  LUX  !-UI 

men             rvta.                            3:  cl_jx^x  iX^qocdc  cxh-  aco 

21  QC      OQQ.DOJ  OZOD-  OO.  Z)l- 

<I  iu      Zh-^cCh-  r£<.~<cL<l  qJLU  aC— t 

CO  h-      'JJcO^QoO  Ol-  I  OLL  QQ  OQ. 


129 


BIBLIOGRAPHY 


1.  Naval  Ship  Research  and  Develooment  Center  Report  No. 
P-433-H-01,  Users  Guiae  NSBDC  Digital  P_£9_3£a.ra  for 
simulating   Submarine  "lotion.    dv   R~onaIcI  W.  Hrcaarcl,  June 

2.  Principles  of  Naval  Architecture,  Revised  Edition,  The 
Society  or~HavaI  IicHitects  ana  Kanne  Engineers. 

3.  Naval  Ship  Research  and  Development  Center  Report  2510, 
Standard  Equations  of  Kotion  for  Sue  sarins  Simulation^  by 
Horton  Gertier  and*  Grant  a.    Hagen,  June  Tvb7. 


4.  Romero,  E.,  Mathematical  Models  and  Computers  Solution 
£2i  tk~  Eguations  or  potion  of  Surface  Shl^s  ancl  Submarines , 
*-Q  2i^  i/e^rees  of  F~s.§iIoj£x.  Masters  Tnesis,  "Naval 
Postgraduate  bcnool7~197Z- 


5«  Fundamentals  of  Submarine  iildrod_ynamicsx  Motionx  and 
Control,  NavshTps  "0"51T-:j'J3-o0"T0. 

6.  Kirk,  D.E.,  Optimal  Control  Theory,  Prentice-Hall  Inc., 
1970. 

7.  Thaler,  G.J.  and  Huang,  J.  ,  Steady-State  Decoupling  and 
Design  of  linear  ^ultiyariable  "systems,  Research  Report, 
University  or  Santa  Clara,  1b>74. 


130 


INITIAL  DISTRIBUTION 

No.  Copies 

1.   Defense  Documentation  Center  2 

Cameron  Station 
Alexandria,  Virginia  22314 


2.   Library,  Code  0212 

Naval  Postgraduate  School 


Monterey,  California  93940 


Department  Chairman,  Code  52 
Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93940 


4,   Professor  G,  Thaler,  Code  52  Tr  10 

Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93940 


A.  Gerba,  Code_52Gz 

Department  of  electrical  Engineering 
Naval  Postgraduate  Scnool 
Monterey,  California  93940 


6.   Professor  D.  KirJc,  Code  52Ki 

Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93940 


7.   Naval  Ship  Research  and  Development  Center 
Annapolis  Laboratory 
Attn;  Wait  Blumberg 
Annapolis,  Maryland  21402 


Department  of  the  Navy 

Naval  Ship  Research  and  Development  Center 

Code  1576 

Attn:  Bill  Smith 

Bethesda,  Maryland  20034 


9.   Department  of  the  Navy- 
Naval  Ship  Research  and  Development  Center 
Code  1576 
Attn:  John  Ware 
Bethesda,  Maryland  20034 


10.  LT  Harold  L.  Drurey,  USN 
39  Allen  Place 
Hartford,  Connecticut  06106 


131 


Thesis 

D777     Drurey 

c.l         Automatic  control  of 

submarine  depth,  pitch 

and  trim. 


!hesD777 

Automatic  control  of 


submarine  depth, 


3"SS2£»8 


